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ABSTRACT 


The  focus  of  this  research  is  to  accurately  model  the  tunnel  junction  interconnect 
within  a  multi-junction  photovoltaic  cell.  A  physically  based  2-D  model  was  created  in 
Silvaco  Inc.’s  software  to  model  the  quantum  tunneling  effect  that  is  realized  within  a 
multi-junction  cell.  The  tunnel  junction  interconnect  is  a  critical  factor  in  the  design  of 
multi-junction  photovoltaics  and  the  successful  modeling  of  the  junction  will  lead  to  the 
ability  to  design  more  efficient  solar  cells.  The  quantum  tunneling  effect  is  based  on  the 
non-local  band-to-band  and  trap  assisted  tunneling  probability  described  by  the  Wentzel- 
Kramers-Brillouin  (WKB)  method. 
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EXECUTIVE  SUMMARY 


The  production  of  electrical  power  is  of  growing  concern  in  today’s  world. 
Renewable  energy  has  been  on  the  forefront  of  this  issue  and  continues  to  provide 
promising  results.  Biofuels,  wind  turbines,  hydroelectric,  and  geothermal  energy  are  all 
sources  of  renewable  energy.  However,  solar  energy  has  proven  to  be  a  viable  candidate 
for  many  niche  applications.  Whenever  there  is  a  need  for  a  small,  lightweight,  low  to  no 
maintenance,  completely  independent  source  of  energy,  photovoltaic  solar  cells  routinely 
provide  a  leading  solution. 

Modern  solar  cell  technology  has  been  around  since  the  1950s  and  has  greatly 
matured  over  the  last  half  century.  Solar  power  has  developed  into  a  multi-billion  dollar 
a  year  industry  [1]  providing  electric  power  to  everything  from  spacecraft  to  automobiles. 
Although  photovoltaic  cells  are  a  proven  technology,  there  is  always  room  for 
improvement.  One  major  disadvantage  of  photovoltaic  technology  is  the  power  density 
of  the  cells,  which  depends  on  the  cell  type  and  cell  location. 

One  solution  to  improve  photovoltaic  power  density  is  to  place  cells  of  different 
material  in  tandem.  These  cells  are  connected  with  a  tunnel  junction  in  order  to  minimize 
shadowing  of  the  sub-cells  and  maximize  the  number  of  photons  absorbed.  According  to 
the  DOE,  a  tunnel  junction  interconnect  is  a  Tier  2  Technical  Improvement  Opportunity 
[1].  The  potential  for  such  a  highly  optimized  photovoltaic  cell  depends  on  the  ability  to 
accurately  model  the  quantum  tunneling  phenomenon  in  the  tunnel  junction. 

The  focus  of  this  research  was  to  model  the  quantum  tunneling  phenomenon  in  a 

© 

degenerately  doped  p-n  junction.  This  was  achieved  using  Silvaco  Inc.’s  ATLAS 
software.  Silicon  (Si),  gallium  arsenide  (GaAs)  and  indium  gallium  phosphide  (InGaP) 
tunnel  junctions  were  realized  using  this  software  and  the  affects  of  varying  doping  levels 
in  response  to  a  forward  bias  was  investigated. 

The  validity  of  the  model  was  tested  by  implementing  the  successfully  developed 
tunnel  junction  into  a  multi-junction  solar  cell  with  a  known  I-V  curve.  The  multi¬ 
junction  cell  model  was  developed  by  Panayiotis  Michalopoulos  [2]  and  the  model  was 


XV 


compared  to  the  experimental  data  found  by  the  Central  Research  Laboratory,  Japan 
Energy  Corporation  [3].  The  simulated  power  density  was  ealculated  to  be  within  3.8% 
differenee  from  the  experimental  results. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Photovoltaic  solar  cells  are  an  exeellent  source  of  renewable  energy  in  countless 
applications.  From  Mars  rovers  to  orbiting  satellites,  the  need  for  a  lightweight, 
renewable,  low  maintenance  energy  souree  is  neeessary  to  extend  equipment  lifetimes 
and  remain  at  the  leading  edge  of  technology  and  exploration.  UAVs,  sensor  nodes  and 
entire  command  posts  are  all  using  solar  energy  to  power  their  systems  to  remain  remote, 
untethered  tools  of  our  military.  Even  the  modern  war  fighter  is  beeoming  more  and 
more  reliant  on  gadget  technology  including  night  vision,  satellite  radios  and  GPS  units, 
whieh  are  routinely  issued  to  provide  an  advantage  over  an  adversary.  Considering  these 
and  many  other  photovoltaie  applications,  efficiency  and  size  are  paramount  in  researeh 
and  development. 

Currently  single  cell  photovoltaic  energy  sources  are  restrieted  to  a  small 
speetrum  of  emitted  energy  depending  on  the  semiconductor  used.  Operating  at  between 
25  to  30%  efficiency,  the  need  for  more  power  only  inereases  array  weight  by  requiring 
more  eells  to  produee  a  specified  output  power.  To  increase  efficieney,  multi-junction 
solar  cells  have  been  developed  in  order  to  absorb  a  wider  array  of  the  solar  spectrum. 
By  stacking  cells  of  eompeting  semiconductor  material,  efficiency  can  theoretically  be 
inereased  to  upwards  of  60%,  whieh  would  lead  to  smaller  arrays  with  increased  power. 

The  one  flaw  to  simply  stacking  cells  one  on  another  is  the  parasitic  junction 
developed  at  the  junetion  between  the  solar  cells.  This  problem  is  solved  by 
manufacturing  a  tunnel  junction  between  each  cell  in  order  to  ensure  that  the  current  flow 
in  the  eell  and  at  the  junction  is  in  the  same  direction. 

B,  OBJECTIVES  AND  APPROACH 

The  purpose  of  this  thesis  is  to  aecurately  model  the  quantum  tunneling 

phenomenon  within  a  degenerately  doped  p-n  junction.  The  junction  will  be  modeled 

using  Silvaco  Data  System  Ine.’s  ATLAS®  software.  ATLAS®  is  a  multi-physics 

semiconductor  simulation  tool  which  provides  the  user  with  the  ability  to  design  and 

1 


evaluate  ehanges  to  a  produet  before  fabrieation.  Onee  designed,  the  developed  junetion 
will  be  applied  to  a  multi-junetion  photovoltaie  eell  in  order  to  inerease  the  effieieney 
while  deereasing  the  development  time  and  money  spent. 

C.  RELATED  WORK 

This  researeh  is  related  to  a  Naval  Postgraduate  Sehool  thesis  eompleted  in  Mareh 
2002  by  Panayiotis  Miehalopoulos  titled  A  Novel  Approach  for  the  Development  and 
Optimization  of  State-of-the-art  Photovoltaic  Devices  Using  Silvaco.  Although  this 
thesis  was  very  sueeessful  in  optimizing  a  multi-junetion  photovoltaie  eell,  the  forward 
bias  I-V  eharaeteristie  of  the  tunnel  junetion  was  not  modeled.  The  ability  to  model  this 
detail  will  further  press  the  optimization  of  photovoltaie  teehnology. 

D,  ORGANIZATION 

The  following  two  ehapters  are  a  diseussion  of  semieonduetor  physios  and  how 
they  apply  to  a  standard  p-n  junetion.  Chapter  IV  extrapolates  this  information  and 
applies  it  to  a  degenerately  doped  tunnel  junetion.  The  simulation  software  along  with 
the  tunnel  junetion  simulation  oode  is  disoussed  in  Chapter  V.  The  results  of  the  tunnel 
junetion  model  are  presented  in  Chapter  VI  followed  by  the  suooessfully  implemented 
dual-junotion  solar  eell  in  Chapter  VII.  Finally,  the  reoommendations  and  oonolusions 
are  presented  in  Chapter  VIII. 
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II.  SEMICONDUCTOR  PHYSICS 


The  basic  material  properties  introduced  in  this  section  are  essential  to  the 
understanding  of  the  following  research.  Although  there  are  more  detailed  works  on  this 
subject,  this  discussion  will  be  limited  to  the  essential  properties  related  to  silicon  (Si), 
germanium  (Ge),  gallium  arsenide  (GaAs)  and  indium  gallium  phosphide  (InGaP); 
semiconductors  typical  to  photovoltaic  applications. 

A,  MATERIAL  SCIENCE  AND  BOHR’S  MODEL 

One  common  way  to  categorize  a  material  is  by  its  electrical  properties. 
Depending  on  the  level  of  an  element’s  resistivity,  it  can  be  categorized  as  an  insulator, 
conductor,  or  semiconductor.  The  ability  of  a  material  to  resist  or  conduct  electricity  is 
dependent  on  many  factors,  including  lattice  structure,  free  electrons,  energy  bandgap, 
and  temperature.  Some  materials  have  very  discrete  electrical  properties  which  define 
them  as  either  an  insulator  or  conductor.  However,  other  materials  such  as  silicon  (Si) 
and  gallium  arsenide  (GaAs)  can  act  as  either  an  insulator  or  conductor  and  are  therefore 
considered  semiconductors.  Figure  1  shows  the  typical  range  of  conductivities  for 
insulators,  conductors,  and  semiconductors. 
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Figure  1.  Resistivity  for  various  material  types,  from  [4]. 


Bohr’s  model  can  be  derived  as  a  first  order  approximation  of  quantum  mechanics 
yet  its  simplicity  is  valuable  to  the  understanding  of  the  electrical  characteristics  of 
materials  on  a  large  scale.  According  to  Niels  Bohr,  atoms  are  comprised  of  three 
subatomic  particles:  a  negative  charged  electron,  a  positive  charged  proton,  and  a  neutron 
with  no  charge.  Compared  to  the  electron,  the  high  mass  proton  and  neutron  form  what  is 
known  as  the  nucleus,  which  sits  at  the  center  of  the  model.  The  electron,  which  is  about 
1,836  times  smaller  than  the  proton,  revolves  around  the  nucleus  at  a  specific  energy 
level  known  as  an  orbit  or  shell.  The  three  subatomic  particles  are  held  together  by  the 
electrostatic  force  between  the  electron  and  proton  similar  to  the  way  planets  are  held  in 
orbit  by  gravity.  Each  shell  must  contain  a  set  number  of  electrons  before  the  outermost 
shell,  known  as  the  valence  shell,  can  be  occupied. 

Modeling  the  simplest  of  atomic  structures,  the  hydrogen  atom,  Bohr  assumes  that 
an  electron  in  orbit  around  the  nucleus  can  take  on  only  certain  levels  of  angular 
momentum  {nh  ).  This  quantization  of  the  electron’s  angular  momentum  is  then  related 

to  the  electron’s  energy  level  by  Equation  2.1  below,  where  m^is  the  mass  of  a  free 
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electron,  q  is  the  electron  charge,  the  permittivity  of  free  space,  n  is  the  electron 


orbit  and  h  is  the  reduced  Plank’s  constant 


V  In 
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Eh=- 


m,q 


13.6 


2i^47rsQn'h) 


eV 


Eqn.  2.1 


As  explained  by  Equation  2.1  and  Eigure  2,  it  is  clear  that  as  the  orbits  increase, 
the  energy  at  that  particular  orbit  also  increases  by  a  factor  of  the  orbit  squared.  It  also 
follows  that  for  an  electron  to  move  from  an  outer  orbit  to  an  inner  orbit  {E^-E^ ),  it 

must  give  off  some  of  its  energy,  which  it  does  in  the  form  of  light  at  a  discrete 
wavelength  (Eqn.  2.2). 


Eigure  2.  Bohr’s  model  of  allowable  electron  orbits  in  a  hydrogen  atom,  from  [5]. 
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Eqn.  2.2 


5 


As  I  will  discuss  later,  this  requirement  of  an  eleetron  to  lose  energy  in  order  to 
move  from  an  outer  shell  to  an  inner  shell  is  the  same  principle  a  solar  cell  uses  to  absorb 
light  and  generate  power. 

B.  CRYSTAL  STRUCTURE 

The  spatial  arrangement  of  the  atoms  within  a  physical  structure  greatly  affects 
the  preeise  properties  of  a  given  material.  As  Figure  3  illustrates,  a  solid  ean  be  elassilfed 
into  three  separate  degrees  of  atomie  order:  amorphous,  polyerystalline  and  crystalline. 
Amorphous  structures  are  arranged  in  such  a  way  that  any  one  section  within  the 
structure  is  different  from  any  other  section.  Crystalline  structures  are  arranged  such  that 
the  atoms  are  arranged  in  the  same  manner  throughout  the  solid.  Polyerystalline 
structures  area  a  eombination  of  both  amorphous  and  crystalline.  Although  it  is  possible 
for  a  semiconductor  to  be  arranged  in  any  of  the  three  atomic  order  classifications, 
polyerystalline  solids  are  most  prominent  in  today’s  market  due  to  cost,  ease  of 
manufaeturing  and  adequate  reliability  in  photovoltaie  applications. 


dancing  y*  ^  '  i 

bonds  ^  <*** J  j 

1  1  1 /'■.v 

1  1  1 
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1  1  1  1  1  1 
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long-range  order 

Completely  ordered 
in  segments 

Entire  solid  is  made  up  of 
atoms  in  an  orderly  array 

Figure  3.  Degree  of  atomic  order  classification  within  solids:  (a)  amorphous  (b) 

polyerystalline  (e)  crystalline,  from  [5]. 
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The  atomic  structure  of  a  semiconductor  can  be  further  described  by  its  lattice 
structure  within  a  unit  cell.  The  typical  way  to  describe  an  atomic  structure  is  to  use  the 
Brevais  lattice  system,  which  describes  the  symmetry  of  a  crystal  by  the  three- 
dimensional  (3-D)  geometric  arrangement  of  the  atoms  in  it.  Although  there  are  seven 
separate  categories  within  the  Brevais  system  with  fourteen  total  arrangements,  the 
simplest  of  the  arrangements  can  be  divided  into  three  categories:  the  simple  cubic,  body- 
centered  cubic  and  face-centered  cubic  (Figure  4). 


Figure  4.  Brevais  lattice  cubics,  from  [6]. 


Taking  a  closer  look  at  the  lattice  structure  of  an  element  with  four  valence 
electrons,  such  as  carbon,  silicon  and  germanium,  there  are  small  differences  from  the 
lattice  structures  described  above  and  those  of  the  group  IV  elements  (those  with  four 
valence  electrons).  Figure  5  is  a  diamond  lattice  in  which  two  face  centered  cubic  lattices 
are  interpenetrating.  The  blue  and  green  spheres  represent  the  first  interpenetrating  face 
centered  cubic  and  the  light  blue  sphere  represents  the  second  interpenetrating  face 
centered  cubic.  The  second  face  centered  cubic  is  diagonally  shifted  Vi  of  the  unit 
distance.  All  element  IV  semiconductors  are  described  as  having  cubic-diamond  lattice 
and  are  so  named  because  each  falls  under  the  carbon  element,  which  is  the  basic 
building  block  of  diamond.  Within  this  lattice  structure  each  atom  is  equidistant  from  its 
four  nearest  neighbors.  Silicon’s  and  germanium’s  lattice  constants  are  5.437  A  and 
5.660  A  respectively. 
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Figure  5.  Cubic-Diamond  lattice  structure,  from  [6]. 


C.  CARRIERS 

Carriers  within  a  material  provide  an  avenue  for  the  transport  of  an  electric  charge 
and  therefore  electric  current.  Similar  to  the  Bohr  model  of  the  hydrogen  atom,  the 
electrons  within  a  Si  atom  only  occur  at  certain  energy  levels.  However,  since  the  silicon 
atom  is  a  multi-electron  atom,  the  associated  shell  model  is  slightly  more  complex.  The 
innermost  filled  shells  along  with  the  nucleus  are  often  referred  to  as  the  core  of  the  atom 
because  they  are  tightly  bound  by  the  electrostatic  force  between  these  shells  and  the 
nucleus.  The  remaining  electrons  of  the  atom  that  fill  the  outer  shell  have  a  weak  bond  to 
the  core,  and  are  referred  to  as  the  valence  electrons  due  to  their  strong  participation  in 
chemical  reactions. 
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Silicon  has  fourteen  electrons,  ten  of  which  occupy  the  two  inner  most  shells 
(atom  eore)  leaving  four  valenee  electrons.  When  the  lattiee  is  formed  between  the 
different  atoms,  there  is  a  finite  probability  that  at  any  temperature  above  zero  Kelvin 
(OK),  an  eleetron  from  the  valenee  band  will  be  knoeked  free  and  an  eleetron  delicieney 
known  as  a  hole  will  be  left  behind.  The  hole  left  behind  retains  a  positive  charge  equal 
to  the  absolute  value  of  the  electron  charge.  The  free  electron  along  with  the  hole  is 
known  as  an  electron  hole  pair  (EHP)  and  can  be  described  by  the  phenomenon  of 
ionization  or  generation.  When  a  free  electron  and  hole  eombine  they  disappear,  a 
proeess  which  is  known  as  recombination.  The  rate  of  reeombination  is  proportional  to 
the  number  of  free  electrons  and  existing  holes.  When  in  thermal  equilibrium  the 
generation  and  reeombination  keep  the  number  of  EHPs  eonstant.  Germanium  reacts  in 
much  the  same  way  although  it  has  thirty-two  total  eleetrons,  twenty-eight  of  whieh 
eonstitute  the  atom  core,  leaving  four  valenee  electrons.  Recombination  and  generation 
will  be  discussed  in  more  detail  later. 

An  intrinsie  material  has  no  impurities  (doping)  in  it  and  will  not  produee  EHPs  at 
OK.  However,  the  EHPs  will  inerease  logarithmically  as  a  function  of  temperature 
(Eigure  6)  and  therefore  conductivity  will  also  increase. 
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Figure  6.  Intrinsic  carrier  density  vs.  Temperature,  from  [7]. 


D.  ENERGY  BANDS 

The  energy  band  model  is  a  useful  way  to  visualize  the  differences  between  an 
insulator,  conductor  and  semiconductor  (Figure  7).  When  looking  at  a  material  rather 
than  a  free  atom,  the  electrons  within  the  material  form  energy  bands  known  as  a  valence 
band  and  conduction  band  at  a  specific  energy  state.  The  difference  between  the  valence 
and  conduction  band  is  known  as  the  bandgap.  The  conductivity  of  a  material  is 
dependent  on  the  number  of  free  electrons  in  the  conduction  band. 

An  insulator  has  a  large  bandgap  that  an  electron  must  overcome  in  order  for  it  to 
flow  in  the  conduction  band.  Conversely,  a  conductor  is  described  as  material  with  the 
valence  and  conduction  bands  overlapping  in  which  case  the  free  electrons  do  not  require 
any  energy  to  traverse  an  energy  barrier.  As  expected,  a  semiconductor  has  a  bandgap 
that  is  not  overlapping  yet  the  energy  difference  between  the  valence  and  conduction 
band  is  minimal. 
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Figure  7.  Energy  band  diagrams,  from  [8] 


Within  the  energy  band  there  is  a  distinetion  between  a  direet  bandgap  and  an 
indireet  bandgap  based  on  the  relative  momentum  between  the  valance  band  hole  and  the 
conduction  band  electron.  A  direct  bandgap  semiconductor  is  one  in  which  the  maximum 
energy  level  of  the  valence  band  and  the  minimum  energy  of  the  conduction  band  occur 
at  the  same  momentum.  Examples  of  direct  bandgap  semiconductors  are  GaAs  and  InP 
which  are  sometimes  referred  to  as  radiant  recombination  semiconductors  because  of  the 
band-to-band  recombination  and  release  of  energy  in  the  form  of  photons  or  radiation 
wavelengths  of  light. 

Indirect  bandgap  semiconductors  have  maximum  and  minimum  energy  band 
levels  at  different  momentums.  Examples  of  indirect  bandgap  semiconductors  are  Si,  Ge 
and  GaP.  The  energy  released  during  indirect  recombination  cannot  be  direct  band-to- 
band  and  requires  the  assistance  of  a  third  party  or  defect  in  the  material.  This  energy  is 
given  off  in  the  form  of  phonons  or  quasi-particles  of  quantized  sound  waves.  The 
energy  band  structure  for  Ge,  Si  and  GaAs  are  compared  in  Eigure  8.  This  relationship 
will  become  more  clear  when  discussing  recombination  and  generation. 
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Figure  8.  Energy  band  strueture  vs.  wave  veetor  for  Ge,  Si  and  GaAs,  from  [7]. 


E,  INTRINSIC  SEMICONDUCTORS  AND  DOPING 

An  intrinsie  semieonduetor  is  one  with  no  impurities  in  it.  Therefore,  it  exhibits 
properties  inherent  to  the  pure  semieonduetor  sample  whieh  ean  be  used  to  identify  the 
material.  Under  equilibrium  eonditions,  an  intrinsie  semieonduetor  has  an  equal  number 
of  holes  (p)  and  eleetrons  (n)  whieh  is  ultimately  equal  to  the  intrinsie  earrier 
eoneentration  (n. ).  The  number  of  p  and  n  are  equal  in  an  intrinsie  semieonduetor 
beeause  earriers  are  only  ereated  in  pairs  in  a  highly  pure  material.  If  an  intrinsie 
semieonduetor  bond  is  broken,  a  free  eleetron  and  a  hole  are  ereated  simultaneously. 
Similarly,  if  an  eleetron  is  exeited  from  the  valenee  band  to  the  eonduetion  band  then  a 
hole  will  simultaneously  be  ereated  in  the  valence  band.  A  list  of  common 
semiconductor  materials  and  their  intrinsic  carrier  concentrations  at  room  temperature  are 
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given  in  Table  1.  Although  these  eoncentrations  may  seem  high,  they  are  only  a  small 
fraetion  of  the  total  number  of  bonds  in  the  respeetive  materials. 


Table  1 .  Intrinsie  earner  eoneentration  at  room  temperature  (300K). 


Carrier  Concentration  ( n^ )  per  crn 

Material 

2x10*^ 

GaAs 

1x10*'’ 

Si 

2x10'" 

Ge 

Adding  impurity  atoms  to  a  material  changes  its  electrical  properties,  a  process 
known  as  doping.  Doping  is  used  to  control  the  number  of  electrons  or  holes  in  a 
material  by  the  addition  of  a  controlled  amount  of  a  specific  impurity  atom.  The  doping 
of  intrinsic  material  results  in  either  a  p-type  or  an  n-type  material,  meaning  it  has  an 
excess  number  of  holes  or  electrons  respectively.  With  respect  to  IV  elements,  acceptors 
(III  elements)  are  added  to  produce  a  p-type  material  and  donors  (V  elements)  are  added 
to  produce  n-type  materials.  This  is  due  to  the  number  of  valence  electrons  in  the  dopant 
atom.  From  Section  II.A  we  know  that  IV  elements  have  four  valence  electrons.  If  an 
element  with  three  valence  electrons  is  added  to  the  diamond  lattice  the  three  valence 
electrons  will  not  be  able  to  complete  all  of  the  semiconductor  bonds  with  the  element  IV 
atoms.  Therefore,  the  excess  electron  will  continue  to  be  accepted  by  the  dopant  atom, 
and  an  excess  hole  will  be  allowed  to  wander  throughout  the  structure  without  the 
creation  of  a  free  electron.  Conversely,  if  a  donor  is  added  to  the  intrinsic  element  IV 
material  then  all  the  semiconductor  bonds  will  be  made,  leaving  an  excess  electron  to 
roam  throughout  the  lattice.  Only  a  small  amount  of  dopant,  on  the  order  of  one  atom  for 
every  10^'’  atoms,  is  required  to  significantly  affect  the  electrical  properties  of  a  material. 
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F.  FERMI  LEVEL 

The  Fermi  level  describes  the  top  of  the  collection  of  electron  energy  levels  at  OK 
and  is  based  on  the  Fermi-Dirac  probability  equation  (Eqn.  2.3).  The  Fermi-Dirac 
probability  equation  is  a  particle  statistic  function  that  determines  the  distribution  of 
fermions  (particles  with  half-integer  spin)  at  a  particular  energy  state  at  absolute  zero. 


f{E)  = 


1 

^(E-E,  )lkT 


Eqn.  2.3 


At  OK  all  the  available  energy  states  below  Ep  will  be  filled  by  fermions  of  one 

and  only  one  particle  which  are  constrained  by  the  Pauli  Exclusion  Principle.  At  low 
temperatures,  any  energy  below  the  Eermi  level  will  have  a  probability  of  1  and  any 
energy  above  the  Eermi  level  will  have  a  probability  of  essentially  0.  As  the  temperature 
rises,  EHPs  are  created  and  electrons  of  higher  energy  than  the  Eermi  energy  emerge  in 
the  conduction  band.  Eigure  9  illustrates  this  distribution  for  increased  temperatures  in  an 
intrinsic  material. 
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Eigure  9.  Eermi  distribution  function,  from  [5]. 


The  Eermi  level  can  be  manipulated  by  doping  the  material  either  with  an 
acceptor  or  donor,  resulting  in  a  p-type  or  n-type  semiconductor  respectively.  By  raising 
or  lowering  the  Eermi  level,  one  can  specify  the  required  energy  for  an  electron  to 
transition  from  the  valance  band  to  the  conduction  band.  As  in  Eigure  10,  an  intrinsic 
material  has  a  Eermi  level  at  the  center  of  the  band  gap.  An  n-type  material  shifts  the 
Eermi  level  closer  to  the  conduction  band  producing  shallow  donors.  Similarly,  a  p-type 
material  shifts  the  Eermi  level  closer  to  the  valence  band  producing  shallow  acceptors. 
With  the  Eermi  level  closer  to  the  valence  and  conduction  bands,  only  a  small  amount  of 
energy  is  required  to  ionize  a  hole  or  electron  into  their  respective  bands. 
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Figure  10.  Schematic  band  diagram,  density  of  states,  Fermi-Dirac  distribution  and  the 
carrier  concentrations  for  intrinsic,  n-type  and  p-type  semiconductors  at  thermal 

equilibrium,  from  [7]. 
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When  the  material  is  not  under  equilibrium  eonditions  the  Fermi  level  is 
influeneed  ereating  Quasi-Fermi  energy  levels  and£'^^)  deseribed  by  the 

distributions  in  Equation  2.4  and  Equation  2.5.  It  is  easy  to  see  from  these  two  equations 
that  at  equilibrium  =  E^^  . 


UE)  = 


1 


Eqn.  2.4 


f,iE)  = 


1 

J^E-E,^]lkT 


Eqn.  2.5 


As  I  will  diseuss  later  with  respeet  to  the  tunneling  effeet,  it  is  also  possible  to 
degenerately  dope  the  material  in  order  to  push  the  Eermi  level  into  the  valenee  or 
eonduetion  band.  However,  it  is  important  to  note  that  a  heavily  doped  material  may 
exhibit  unintended  eleetrieal  properties  depending  on  the  distribution  of  the  dopant. 
Beeause  the  distribution  of  the  dopant  may  not  be  homogenous  throughout,  different 
energy  fluetuations  may  be  found  at  different  positions  within  a  strongly  eompensated 
semieonductor.  This  is  due  to  the  mismateh  of  the  lattiees  between  the  two  materials, 
whieh  ean  greatly  affeet  the  eonduetivity  of  a  material. 

G.  MOBILITY 

Mobility  (p)  is  the  ratio  of  average  speed  (drift  speed)  to  the  applied  field  of  a 
material  and  is  analogous  to  permittivity  ).  Ideally  (in  a  vacuum)  an  electron  inside  of 

an  electric  field  will  continuously  accelerate.  However,  in  a  material  the  electron  will 
accelerate  only  until  it  collides  with  an  atom  in  the  lattice.  As  expected,  the  higher 
doping  in  a  material  will  result  in  a  higher  probability  of  a  collision  with  atoms  in  the 
lattice  and  therefore  mobility  will  decrease.  Eigure  11  illustrates  this  point  for  some 
common  semiconductors. 
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Figure  1 1 .  Mobility  vs.  doping  level,  from  [7]. 


H,  CURRENT  DENSITY 

The  total  particle  current  in  a  material  is  equal  to  the  total  current  density  (J)  of 
both  the  holes  ( Jp )  and  electrons  ( )  which  results  from  the  drift  and  diffusion  current 
densities.  Diffusion  current  (I^,)  is  the  result  of  electrons  or  holes  in  a  concentrated  area 
spreading  themselves  evenly  though  a  material  (Figure  12)  in  order  to  achieve  a  zero 
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concentration  gradient  (V«  =  0  for  holes,  V/7  =  0  for  eleetrons)  throughout.  Diffusion 
eurrent  density  ( for  holes,  Jj^|j;j.pfor  eleetrons)  is  analogous  to  the  majority  earrier 

eoncentration  (doping)  of  a  material  and  is  deseribed  by  Equation  2.6  and  Equation  2.7 
where  Dp  and  are  the  eonstants  of  proportionality. 


Id 

- ► 


Jp|d.ff  =  -qD^'^P  Eqn.  2.6 

JNidiff  =  Eqn.  2.7 


Drift  eurrent  ( Ig )  is  the  result  of  the  introduetion  of  an  eleetrie  or  magnetie  field 

and  is  analogous  to  the  minority  earrier  eoneentration  and  is  proportional  to  the  intensity 
of  the  field  (rf).  Eigure  13  illustrates  drift  eurrent  in  the  presenee  of  an  eleetrie  field  and 
Equation  2.8  and  Equation  2.9  deseribe  the  mathematical  model  of  drift  eurrent  as  a 
function  of  mobility  (p)  and  applied  eleetrie  field. 
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Figure  13.  Drift  current  in  the  presence  of  an  electric  field,  from  [2]. 

Eqn.  2.8 
Eqn.  2.9 


Jpldrift 
JNldrift  = 


Is 

► 


Both  the  electron  ( //^ )  and  hole  ( //^ )  current  densities  depend  on  their  associated 

diffusion  and  drift  currents  (Eqn.  2.10  and  Eqn.  2.11)  which  ultimately  lead  to  the  total 
current  density  (Eqn.  2.12). 


■^P  "^P  drift  "^Pdiff  qD^Vp 

Eqn.  2.10 

■^N  “  -^N  drift  "*"'^Ndiff  “ 

Eqn.  2.1 1 

J  —  +  Jp 

Eqn.  2.12 

I.  GENERATION  AND  RECOMBINATION 

Generation  is  the  process  whereby  electrons  and  holes  (carriers)  are  created. 
Conversely,  recombination  is  the  process  whereby  electrons  and  holes  are  annihilated. 
Both  generation  and  recombination  are  order-restoring  mechanisms  to  return  the  system 
to  its  equilibrium  state  (p=n=  n. )  after  the  carrier  excess  or  deficit  in  the  semiconductor  is 
stabilized  or  eliminated. 
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There  are  three  important  recombination-generation  methods  as  shown  in  Figure 
14.  Band-to-band  recombination  (a)  is  the  direct  annihilation  of  a  valance  band  hole  and 
a  conduction  band  electron  when  the  hole  and  electron  stray  into  the  same  vicinity  within 
the  semiconductor  lattice.  The  loss  of  energy  experienced  by  the  electron  must  release  an 
energy  quantum  equal  to  the  total  energy  lost.  This  loss  is  typically  realized  in  the  form 
of  light  as  described  in  Equation  2.2. 

The  R-G  (recombination-generation)  center  recombination  is  caused  by  defects  or 
special  impurity  atoms  within  the  lattice  and  takes  place  only  in  special  locations  within 
the  semiconductor  known  as  R-G  centers.  The  R-G  center  concentration  is  typically  very 
low  compared  to  the  donor  and  acceptor  concentration  yet  leads  to  the  introduction  of 
various  energy  levels  near  the  center  of  the  band  gap  known  as  the  trap  energy  {Ej.). 

The  R-G  center  recombination  is  a  two-step  process  and  can  take  place  in  one  of  two 
ways  as  represented  in  Figure  14(b).  The  first  method  is  when  a  carrier  (hole  or  electron) 
is  trapped  at  the  impurity  location  until  the  opposite  carrier  is  attracted  to  the  first  and 
recombination  annihilates  them  both.  The  second  is  when  an  electron  falls  into  the 
impurity  location  where  it  will  initially  lose  energy  equal  to  Once  a  hole 

comes  into  vicinity,  the  electron  is  attracted  into  the  valence  band  after  which  it  will  lose 
the  rest  of  its  energy  equal  to  \Ej.  -  Ey  |  before  it  recombines  with  the  hole  and  disappears. 

This  type  of  recombination  will  cause  quantum  lattice  vibrations  which  will  release 
thermal  energy  and  ultimately  increase  the  temperature  of  the  semiconductor. 

The  third  process  is  called  Auger  recombination  (c),  which  is  also  a  band-to-band 
recombination.  However,  this  process  occurs  simultaneously  with  a  collision  between 
two  like  carriers.  The  energy  transferred  between  the  colliding  carriers  forces  one  carrier 
into  the  alternate  energy  band  while  the  second  carrier  remains  in  the  original  band.  The 
second  carrier  becomes  highly  energetic  and  then  returns  to  its  original  energy  level  in 
small  heat  producing  steps  as  it  collides  with  the  semiconductor  lattice.  This  unique 
stepwise  energy  loss  or  “staircase”  is  defined  as  thermalization  and  also  leads  to  an 
increased  semiconductor  temperature. 
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Figure  14.  Energy  band  visualization  of  recombination  and  generation  processes, 

from  [5]. 


The  generation  process  is  simply  the  reverse  of  the  recombination  process.  For 
band-to-band  generation  (d),  thermal  energy  or  light  energy  >^g)  absorbed  by  a 

semiconductor  will  excite  electrons  from  the  valence  band  into  the  conduction  band 
creating  an  EHP.  Similarly,  for  R-G  center  generation  (e),  thermal  energy  absorbed  by 
the  semiconductor  will  excite  the  electron  enough  to  achieve  the  trap  energy  at  the 
impurity  location  and  then  again  to  the  conduction  band.  Carrier  generation  via  impact 
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ionization  (f)  is  analogous  to  Auger  recombination,  because  the  creation  of  EHPs  is  the 
result  of  energy  released  when  a  highly  energetic  carrier  collides  with  the  semiconductor 
lattice. 


J.  SUMMARY 

The  basic  understanding  of  semiconductor  physics  leads  to  a  more  advanced 
discussion  of  a  p-n  junction,  which  is  discussed  in  the  next  chapter. 
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III.  P-N  JUNCTIONS 


A,  INTRODUCTION  TO  JUNCTIONS 

The  successful  implementation  of  a  P-N  junction  is  critical  in  modern  electronics 
and  semiconductor  devices,  including  bipolar  junction  and  metal  on  oxide  transistors 
(BIT  and  MOSFET),  rectifier  and  Zener  diodes,  and  photovoltaics  to  name  a  few. 
Through  the  doping  process,  a  P-N  junction  is  formed  to  effectively  control  the  electrical 
characteristics  of  a  material. 


B.  BASIC  STRUCTURE 

A  p-n  junction  is  formed  by  a  single  crystal  semiconductor  where  one  side  is 
doped  with  acceptors  (p)  and  the  other  side  is  doped  with  donors  (n).  At  the  location  of 
contact  the  metallurgical  junction  is  formed  as  shown  in  Figure  15.  Initially  there  is  a 
very  large  majority  carrier  density  gradient  at  the  metallurgical  junction.  Equilibrium 
will  be  sought  as  the  majority  carriers  in  the  p-type  material  (holes)  will  begin  to  diffuse 
into  the  n-type  material  while  the  majority  carriers  from  the  n-type  material  (electrons) 
will  begin  to  diffuse  into  the  p-type  material. 


\ 


Metallurgical 

junction 


Hole  diffusion 

Electron  diffusion 

x  =  0 

(a)  (b) 

Figure  15.  (a)  Simplified  P-N  junction;  (b)  doping  profde,  from  [9]. 
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The  diffusing  carriers  from  both  regions  will  combine  with  the  minority  carriers 
on  the  opposite  side  and  disappear.  Assuming  that  there  is  no  outside  electric  field 
connected  to  the  junction,  this  diffusion  can  not  continue  indefinitely.  The  holes 
diffusing  from  the  p-region  will  leave  behind  negatively  charged  acceptor  atoms  and  the 
electrons  diffusing  from  the  n-region  will  leave  behind  positively  charged  donor  atoms. 
The  net  charges  in  the  p-  and  n-regions  will  induce  an  electric  field  in  the  metallurgical 
junction,  which  will  form  what  is  known  as  a  space  charge  region  or  depletion  region  as 
shown  in  Figure  16. 


Na  negative 
charge 


Nd  positive 
charge 


■< -  Space  charge  region  - ► 


“Diffusion 
focus"  on 
holes 


E-field 


E-field 
•«-  force  on 
holes 


E-field 
force  on  -► 
electrons 


Figure  16.  P-N  junction  depletion  region  formed  by  diffusion,  from  [9]. 


When  in  thermal  equilibrium,  the  Fermi  levels  on  both  sides  of  the  junction  are 
equal  (Figure  17).  The  potential  barrier  known  as  a  built  in  voltage  (F^.)  is  responsible 
for  maintaining  equilibrium  between  the  majority  carrier  electrons  in  the  n-region  and  the 
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minority  carrier  electrons  in  the  p-region.  Similarly  for  the  holes  in  the  valence  band,  the 
built  in  voltage  maintains  equilibrium  between  the  majority  carrier  holes  in  the  p-region 
and  the  minority  carrier  holes  in  the  n-region.  As  expected,  the  built  in  voltage  depends 
on  doping  levels  ( and  A^)  within  the  n-  and  p-regions  (Eqn.  3.1). 
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Figure  17.  P-N  junction  energy  band  diagram  with  Fermi  level. 


7/ 

Vt,  = — In 

q 


Eqn.  3.1 


The  width  of  the  depletion  region  (W)  is  also  related  to  the  relative  doping  level 
of  the  two  regions.  This  characteristic  will  become  more  important  with  the 
implementation  of  the  tunnel  junction  (Section  IV).  The  width  of  the  depletion  region  (in 
cm)  can  be  calculated  by  using  Equation  3.2. 


W  = 


i  ^ 

_  Jj 

Eqn.  3.2 
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Applying  different  bias  eonditions  between  the  two  regions  will  disturb  the  p-n 
junetion  from  equilibrium  and  the  Fermi  levels  will  no  longer  remain  eonstant.  For  a 
standard  doped  junetion,  the  three  conditions  of  interest  are  reverse  bias,  forward  bias  and 
the  breakdown  voltage. 


C.  REVERSE  BIAS 

Applying  a  negative  voltage  (F^)  to  the  p-region  with  respect  to  the  n-region 

results  in  a  reverse  bias  condition  as  shown  in  Figure  18.  The  voltage  applied  attracts  the 
majority  carriers  from  their  respective  regions  and  leaves  minority  carriers  in  its  place. 
This  increases  the  voltage  barrier  and  widens  the  depletion  region.  The  difference 
between  the  two  Fermi  levels  is  equal  to  the  bias  voltage.  Therefore,  the  total  barrier 
voltage  ( )  realized  by  the  system  is  the  sum  of  the  built  in  voltage  and  bias  voltage 
(Eqn.  3.3).  This  barrier  is  easiest  to  visualize  by  the  use  of  an  energy  band  diagram  in 
Figure  19. 
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Figure  18.  Reversed  bias  p-n  junction  schematic,  from  [2]. 
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Eqn.  3.3 
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Eigure  19.  Reversed  bias  p-n  junction  energy  band  diagram. 


The  depletion  region  width  is  given  by  Equation  3.4.  When  in  reverse  bias  the 
probability  that  an  electron  will  gain  enough  energy  to  overcome  the  potential  barrier  is 
very  small  and  the  junction  will  act  like  an  insulator,  not  allowing  current  to  flow.  With 
the  increased  depletion  region,  diffusion  current  {Id)  is  very  small  and  the  total  current  in 
the  system  is  approximately  equal  to  the  drift  current  {!$)  (Eqn.  3.5).  Ideally,  this  current 
would  be  equal  to  zero,  but,  in  reality,  there  is  a  small  amount  of  drift  current  in  the 
system,  as  illustrated  in  Eigure  20. 


W  = 


\2e,(r„+K) 
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Eqn.  3.4 
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I  total  =Is+Id~  Is  (very  small) 


Eqn.  3.5 


Figure  20.  Reversed  bias  p-n  junction  I-V  characteristic. 

D.  FORWARD  BIAS 

Applying  a  positive  voltage  (F^)  to  the  p-region  with  respect  to  the  n-region  is 

considered  forward  bias  (Figure  21).  The  forward  bias  opposes  the  built  in  field 
developed  in  the  depletion  region  and  the  region  becomes  narrower.  When  the  junction 
is  forward  bias,  the  Fermi  level  in  the  p-region  is  increased  and  the  Fermi  level  in  the  n- 
region  is  decreased  to  a  point  where  electrons  can  overcome  the  small  potential  barrier 
(Eqn.  3.6),  as  illustrated  by  Figure  22. 
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Figure  21 .  Forward  bias  p-n  junction  schematic,  from  [2], 


Eqn.3.6 

Ito,ai=h+^D~lD  Eqn.  3.7 
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Figure  22.  Forward  bias  p-n  junction  energy  band  diagram. 


As  the  continuously  supplied  majority  carriers  from  the  n-region  transition  into 
the  p-region,  they  become  minority  carriers  where  they  diffuse  and  eventually  combine 
with  the  majority  holes.  Similarly,  the  continuously  generated  holes  in  the  p-region  will 
be  available  to  combine  with  the  electrons  traversing  the  potential  barrier.  Just  as  before, 
some  of  the  current  being  produced  is  due  to  drift  (Is)  and  some  is  due  to  diffusion  (Id). 
However,  in  the  forward  bias  condition  the  diffusion  current  is  much  greater  than  the  drift 
current  (Eqn.  3.7)  due  to  the  decreased  depletion  region.  Figure  23  shows  the  I-V 
characteristic  of  a  p-n  junction  in  the  forward  bias  condition.  Notice  that  the  current  does 
not  begin  to  flow  until  sufficient  voltage  is  applied  to  overcome  the  potential  barrier. 
This  happens  at  approximately  0.7V  for  Si  and  O.IV  for  Ge. 
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Figure  23.  Forward  bias  p-n  junction  I-V  characteristic. 


E.  BREAKDOWN/ZENER  TUNNELING 

Adding  reverse  bias  beyond  a  certain  threshold  voltage  will  result  in  the 
breakdown  of  the  p-n  junction.  This  breakdown  phenomenon  is  due  to  either  a  Zener 
effect  or  an  avalanche  effect. 

The  Zener  or  tunneling  effect  is  similar  to  the  tunneling  which  I  will  discuss  at 
length  in  the  next  section.  The  difference  here  is  that  the  Zener  effect  takes  place  when  a 
minority  carrier  in  the  p-region  valence  band  is  able  to  penetrate  through  the  energy 
barrier  to  an  empty  state  in  the  n-region  conduction  band  unaffected.  This  can  only  occur 
once  the  valence  band  of  the  p-region  increases  beyond  the  conduction  band  of  the  n- 
region  due  to  the  reverse  bias.  The  exchange  of  minority  carriers  creates  a  large  change 
in  reverse  current  with  a  small  change  in  reverse  bias.  Zener  diodes  are  designed  to 
achieve  this  result  (Figure  24). 
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Figure  24.  Zener  breakdown  p-n  junetion  I-V  eharaeteristie. 


The  avalanche  effect  is  also  a  minority  carrier  effect  in  which  an  electron  from  the 
p-region  with  a  very  large  kinetic  energy  collides  with  atoms  and  breaks  the  covalent 
bonds.  Often  the  EHPs  created  collide  with  other  atoms  to  create  even  more  EHPs, 
creating  an  avalanche-type  effect.  As  more  and  more  minority  carriers  are  created,  they 
are  swept  to  the  opposite  region  which  creates  a  large  reverse  current.  This  runaway 
effect  often  leads  to  excess  heat  from  the  lattice  vibrations  and  results  in  a  damaged  p-n 
junction. 

F.  LOSSES 

Ideally  there  would  be  no  losses  in  the  system  and  the  slope  of  the  I-V  curve 
would  be  perpendicular  to  the  x-axis  when  on.  However,  since  the  semiconductor 
material  has  an  inherent  ohmic  resistance  due  to  atomic  masses  and  impurities,  there  is  a 
slight  ohmic  loss  in  the  junction.  Putting  three  regions  together,  a  typical  p-n  junction  I- 
V  curve  is  realized  (Eigure  25). 
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Figure  25.  p-n  junction  I-V  characteristic. 


G.  SUMMARY 

Controlling  the  doping  levels  in  the  respective  regions  will  shift  the  Fermi  levels 
and  affect  the  ability  of  an  electron  or  hole  to  move  through  the  lattice.  Increasing  the 
doping  levels  to  a  point  where  the  Fermi  levels  are  within  the  energy  bands  of  the  n-  and 
p-  regions  is  essential  to  tunneling  phenomenon  and  will  be  discussed  in  the  following 
Chapter. 
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IV.  TUNNEL  JUNCTION 


A.  TUNNELING  THEORY 

The  tunneling  phenomenon  is  a  majority  earner  effect  in  which  both  the  n-  and  p- 
regions  are  degenerately  doped  (very  heavily  doped  with  impurities)  allowing  an  electron 
to  tunnel  through  a  potential  barrier  that,  under  normal  conditions,  would  not  be  possible 
to  penetrate.  The  tunneling  is  based  on  a  quantum  transition  probability  (7^)  per  unit 
time  and  is  proportional  to  the  exponent  of  the  average  value  of  momentum  ( ^  (o) )  within 
the  tunneling  path  across  the  depletion  region  (Eqn.  4.1).  This  momentum  corresponds  to 
an  incident  carrier  with  a  transverse  momentum  of  zero  and  energy  equal  to  the  Fermi 
energy.  The  tunneling  probability  equation  will  be  further  explained  later  in  this  section. 


7;  =  exp 


-2[  "  A:(x)  dx 

J-x, 


Eqn.  4.1 


B,  PURPOSE 

To  achieve  a  higher  power  density,  multiple  cells  with  different  I-V  curves  and 
decreasing  bandgaps  are  placed  in  tandem  as  illustrated  in  Figure  26.  As  a  result  the 
incident  light  will  be  absorbed  by  the  top  cell  allowing  the  unabsorbed  photons  to  pass 
through  to  the  sub  cells.  This  results  in  the  sum  of  the  voltage  produced  by  the  individual 
cells. 
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Figure  26. 


Tandem  solar  cells. 


Unfortunately,  stacking  the  cells  in  tandem  would  result  in  a  parasitic  junction 
formed  by  the  two  dissimilar  regions  of  the  different  cells  (Figure  27)  and  must  be 
separated  in  order  to  allow  current  flow. 
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Figure  27.  Tandem  solar  cell  with  parasitic  junction. 
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One  option  is  to  mechanically  stack  the  cells  using  an  oxide/conductor  contact  as 
shown  in  Figure  28.  However,  the  added  contact  material  shadows  the  sub  cells  from  the 
unused  photons  being  passed  through  the  top  cell,  limiting  any  efficiency  gained  by 
placing  the  cells  in  tandem. 
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Figure  28.  Tandem  solar  cell  example  with  oxide/conductor  contact  separation,  from  [2]. 

Another  solution  to  this  issue  is  to  use  a  tunnel  junction  to  separate  the  individual 
cells  in  the  tandem  cell.  The  tunnel  junction  will  allow  photons  to  pass  through  as  well  as 
pass  current  with  minimal  voltage  loss  (Figure  29). 
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Figure  29.  Tandem  cell  with  tunnel  junction. 


C.  THE  JUNCTION 

A  tunnel  junction  is  simply  a  tunnel  diode  or  degenerately  doped  p-n  junction 
commonly  used  in  switching  and  oscillator  circuit  applications.  The  difference  between 
the  tunnel  junction  and  the  tunnel  diode  is  that  the  tunnel  junction  is  grown  as  a  single 
crystal  joining  two  solar  cells,  whereas  a  tunnel  diode  operates  as  a  discrete  component 
within  a  circuit.  The  purpose  of  mentioning  the  tunnel  diode  is  for  completeness,  as 
tunnel  diodes  have  been  studied  and  used  in  electronics  since  their  invention  by  Leo 
Esaki  in  1958,  who  received  the  Nobel  Prize  in  Physics  in  1973  for  his  work.  From  this 
point  on,  there  will  be  no  distinction  made  between  the  two  and  I  will  focus  on  the  tunnel 
junction  and  its  application  to  photovoltaics. 

A  tunnel  junction  is  formed  with  degenerately  doped  regions  in  order  to  push  the 
respective  Fermi  levels  into  the  valence  and  conduction  bands  (Figure  30).  As  a  result, 
the  p-region  valence  band  will  contain  holes  and  the  n-region  conduction  band  will 
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contain  electrons  at  thermal  equilibrium.  The  depletion  region  width  is  also  affeeted  by 
the  degenerately  doped  regions  eompressing  it  to  <100A,  whieh  is  essential  to  the 
tunneling  effeet  and  is  signifieantly  smaller  than  a  eonventional  p-n  junetion. 
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Figure  30.  Degenerately  doped  tunnel  junction  energy  band  diagram. 


In  order  for  quantum  tunneling  to  occur,  four  criteria  must  be  met.  The  necessary 
conditions  for  tunneling  are:  (1)  oecupied  energy  states  exist  on  the  side  from  whieh  the 
electron  will  tunnel;  (2)  unoccupied  energy  states  exist  at  the  same  energy  levels  as  those 
electrons  in  the  occupied  energy  states;  (3)  the  tunneling  barrier  width  is  small  and  the 
potential  barrier  height  is  low  enough  that  there  is  a  finite  tunneling  probability;  and  (4) 
the  momentum  is  eonserved  in  the  tunneling  proeess  [2]. 

D.  THERMAL  EQUILIBRIUM 

As  in  the  eonventionally  doped  p-n  junetion,  the  two  Fermi  levels  of  the  p-  and  n- 
regions  are  equal  at  thermal  equilibrium.  Figure  31  illustrates  the  tunneling  junetion 

energy  bands  at  thermal  equilibrium  as  well  as  the  eorresponding  I-V  plot.  The  red  dot 
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on  the  I-V  plot  corresponds  to  the  location  of  interest  on  the  curve.  Since  both  regions 
are  doped  to  the  point  where  the  Fermi  levels  exist  inside  the  valence  and  conduction 
band  and  are  therefore  considered  degenerately  doped. 


Figure  3 1 .  Tunnel  junction  energy  band  and  I-V  characteristic  at  thermal  equilibrium, 

from  [7]. 


E.  REVERSE  BIAS 

Similar  to  the  Zener  effect  discussed  in  the  last  section,  in  reverse  bias  the 
minority  carrier  in  the  p-region  valence  band  is  able  to  tunnel  through  the  potential 
barrier  to  the  conduction  band  of  the  n-region.  Unlike  the  Zener  effect  in  a  standard  p-n 
junction,  the  reverse  bias  tunneling  in  this  junction  occurs  as  soon  as  a  reverse  bias  is 
applied  to  the  junction.  The  reason  for  the  immediate  response  is  the  degenerate  doping 
of  the  regions.  Since  the  Fermi  levels  are  inside  the  respective  valence  and  conduction 
bands  at  equilibrium,  there  are  empty  states  available  for  the  minority  carrier  to  fill  as 
soon  as  reverse  bias  is  applied.  Therefore,  a  large  reverse  current  is  realized  for  a  small 
change  in  voltage.  Figure  32  illustrates  the  minority  carrier  electron  tunneling  to  the 
conduction  band  of  the  n-region  as  well  as  the  location  on  the  I-V  curve  where  reverse 
bias  tunneling  takes  place. 
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Figure  32.  Tunnel  junction  energy  band  and  I-V  characteristic  in  forward  bias 

(reverse  bias),  from  [7]. 


F.  FORWARD  BIAS 

At  forward  bias  (Figure  33),  the  tunnel  junction  has  available  energy  states  in  the 
p-region  corresponding  to  filled  energy  states  in  the  n-region.  As  the  voltage  is  increased 
from  zero  to  ,  the  valence  band  energy  of  the  p-region  and  the  conduction  band  energy 

of  the  n-region  become  closer  in  value. 
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Figure  33 .  Tunnel  junction  energy  band  and  I-V  characteristic  in  forward  bias 

(peak  tunneling),  from  [7]. 


G.  NEGATIVE  RESISTANCE 

From  Fp  to  ,  tunneling  is  still  occurring  but  fewer  and  fewer  empty  states  are 
available.  With  a  small  increase  in  forward  bias  (>  ),  the  number  of  electrons  in  the  n- 

region  directly  opposite  the  empty  states  in  the  p-region  continues  to  decrease  which  will 
result  in  a  decreased  tunneling  current.  The  region  from  to  is  described  as 

negative  resistance  (Figure  34)  and  will  continue  until  the  valence  band  of  the  p-region  is 
equal  to  the  conduction  band  of  the  n-region  at  which  time  the  tunneling  current  will  be 
zero. 
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Figure  34.  Tunnel  junction  energy  band  and  I-V  characteristic  in  forward  bias 

(negative  resistance),  from  [7]. 


H,  THERMAL  CURRENT  FLOW 

Increasing  the  forward  bias  even  further  beyond  will  eventually  result  in  the 

normal  current  flow  of  a  p-n  junction.  As  illustrated  in  Figure  35,  the  energy  band 
diagrams  of  the  p-  and  n-regions  are  reminiscent  of  the  standard  p-n  junction  with  an 
applied  forward  bias.  The  relationship  between  a  standard  p-n  junction  and  a 
degenerately  doped  tunnel  junction  can  be  further  understood  by  comparing  the  two 
curves  in  Figure  36.  As  additional  forward  bias  is  applied,  the  current  flow  of  the 
standard  p-n  junction  and  that  of  the  tunnel  junction  merge. 
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Figure  35 .  Tunnel  junction  energy  band  and  I-V  characteristic  in  forward  bias  (thermal 

current  flow),  from  [7]. 


Figure  36.  Characteristic  I-V  curve  of  a  tunnel  (red)  and  conventional  (blue)  p-n  junction. 
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I,  DIRECT  AND  INDIRECT  TUNNELING 

Just  as  there  is  a  direct  and  indirect  materials  relationship,  so  too  is  there  a  direct 
and  indirect  tunneling  process.  Superimposing  the  E-k  relationship  at  the  classical 
turning  points  onto  the  E-x  relationship  of  the  tunnel  junction  reveals  the  direct  and 
indirect  tunneling  differences.  In  Eigure  37  (a)  and  (b)  the  electrons  in  the  vicinity  of  the 
minimum  of  the  conduction  band  energy-momentum  can  tunnel  to  a  corresponding  value 
of  momentum  in  the  vicinity  of  the  valence  band.  However,  the  distinction  lies  in  the 
electron  and  hole  momentum. 

Eor  direct  band  tunneling  to  occur,  the  momentum  for  the  valence  band  maximum 
and  the  conduction  band  minimum  must  be  the  same.  This  type  of  tunneling  is  illustrated 
in  Eigure  37  (a)  and  takes  place  in  direct  bandgap  semiconductors  such  as  GaAs  and  InP 
as  well  as  indirect  bandgap  semiconductors  where  the  voltage  applied  to  the  junction  is 
large  enough  that  the  valence  band  maximum  (E  point)  is  in  line  with  the  conduction 
band  minimum  ( E  point)  [7]. 

Conversely,  indirect  band  tunneling  occurs  when  the  conduction  band  minimum 
and  the  valence  band  maximum  do  not  occur  at  the  same  momentum,  as  shown  in  Eigure 
37  (b).  As  a  result,  momentum  is  conserved  by  scattering  agents  such  as  photons  or 
impurities  and  is  deemed  photon  assisted  tunneling.  Due  to  the  complex  nature  of  photon 
assisted  tunneling,  it  is  far  less  likely  to  take  place  when  direct  tunneling  is  possible. 
However,  it  is  possible  for  both  direct  and  indirect  tunneling  to  occur  at  the  same  time. 
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Figure  37.  (a)  Direct  band  tunneling  process,  (b)  Indirect  band  tunneling  process, 

from  [7] 


J.  TUNNELING  CURRENT 

For  quantum  tunneling  to  exist,  the  electric  field  in  the  depletion  region  must  be 
on  the  order  of  lO^V/cm.  When  this  is  achieved,  there  is  a  finite  probability  that  an 
electron  will  be  excited  from  the  valence  band  directly  to  the  conduction  band  of  the 
semiconductor. 

Local  band-to-band  tunneling  models  in  Silvaco’s  ATLAS  account  for  the 
electric  field  at  each  node  to  calculate  the  generation  rate  due  to  tunneling.  However,  in 
reality  the  tunneling  is  non-local  and  it  is  therefore  necessary  to  take  the  spatial  profile  of 
the  energy  bands  into  account  as  well  as  the  spatial  separation  of  the  electrons  generated 
in  the  conduction  band  from  the  holes  generated  in  the  valence  band  [9].  To  achieve  this, 
ATLAS  creates  a  series  of  1-D  slices  through  the  device  to  calculate  the  tunneling 
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current.  The  slices  are  generated  as  an  additional  mesh  through  a  series  of  qty  and  qtx 
statements,  which  are  then  plaeed  over  the  top  of  the  original  mesh.  We  will  diseuss  the 
mesh  statements  in  detail  in  Seetion  V. 

With  sufficiently  high  doping  of  both  the  p-  and  n-region  of  the  junction,  any 
electron  within  the  appropriate  energy  range  has  a  finite  probability  of  tunneling  from  the 
valenee  band  to  the  conduction  band.  ATLAS  considers  eaeh  energy  level  within  this 
range  and  ealeulates  the  spatial  start  ( -Xj )  and  end  ( Xj )  point  at  that  level.  The  tunneling 

eurrent  ean  then  be  ealeulated  as  a  funetion  of  the  tunneling  probability  (T  (£')),  the 

effective  mass  (m*),  and  the  potential  barrier  height  {E^)  as  deseribed  by  Equation  4.2, 

where  k  is  Boltzmann’s  eonstant,  q  the  eleetrie  charge,  T  is  temperature,  ti  is  the  redueed 
planks  eonstant  and  E  is  the  eleetron  energy. 


J(E) 


^  ’  [l  +  exp[E^^-E]/kr 


AE 


Eqn.  4.2 


ATEAS  provides  the  tools  to  implement  a  number  of  probability  equations  for 
both  loeal  and  non-local  tunneling  but  none  proved  better  at  aehieving  the  forward  bias 
than  the  Wentzel-Kramers-Brillouin  method  (WKB). 

K.  WENTZEL-KRAMERS-BRILLOUIN  METHOD  (WKB) 

The  tunneling  probability  is  given  by  the  WKB  method  (Eqn.  4.3)  in  whieh  the 
wave  vector  (A:(x))  of  a  carrier  is  evaluated  as  an  exponential  function  and  integrated 
over  the  depletion  region  from  -Xj  to  the  Xj . 


□  exp 


-ij  "  |A:(x)  dx 


Eqn.  4.3 
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Although  there  are  many  ways  to  implement  this  equation  based  on  the  geometry 

of  the  potential  barrier,  the  method  I  will  use  is  based  on  a  triangular  geometry  (Figure 

(6) 

38)  because  this  is  the  method  used  by  the  ATLAS  modeling  software. 


Figure  38.  Triangular  potential  barrier  for  tunneling  model,  from  [7] 


For  the  triangular  potential  barrier  in  Figure  38  the  wave  vector  is  given  as  [7]; 


A:(x) 


2m 


[PE-E)  = 


2m 


-  q/oc 


Eqn.  4.4 


Where  PE  is  the  potential  energy,  E  is  the  incoming  electron  energy.  Eg  is  the  energy 
bandgap,  and  F”  is  the  electric  field. 


Substituting  Equation  4.4  into  Equation  4.3  yields  [7]; 
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The  quantum  tunneling  probability  ean  be  simplified  to: 


4^ 

□  exp - - -  Eqn.  4.6 

3  qf^h 


Where  is  the  bandgap  of  the  material,  m*  is  the  effeetive  mass  defined  in 
Equation  4.7,  q  is  the  eleetrie  eharge,  g’  is  the  eleetrie  field,  and  fi  is  the  redueed 
Plank’s  eonstant. 

m*  =  (-^2)^/2  (-^i)  Eqn.  4.7 


51 


L,  SUMMARY 

The  tunneling  phenomenon  presented  is  based  on  the  probability  that  an  electron 
is  able  to  penetrate  a  potential  barrier  if  there  are  occupied  and  unoccupied  states  at  the 
same  energy  level  on  both  sides  of  the  potential  barrier,  the  distance  between  the 
beginning  and  end  of  the  potential  barrier  is  sufficiently  small,  and  momentum  of  the 
tunneling  electron  is  conserved.  The  following  Chapter  will  implement  the  given 
equations  into  the  Silvaco  software. 


52 


V.  SIMULATION  SOFTWARE 


A.  SILVACO 

Silicon  Valley  Company  (Silvaco)  is  a  leading  vendor  in  technology  computer 
aided  design  (TCAD).  Established  in  1984  and  located  in  Santa  Clara,  California, 
Silvaco  has  developed  a  number  of  exeeptional  CAD  simulation  tools  to  aid  in 
semiconduetor  proeess  and  device  simulation.  Silvaco  has  an  extensive  support  team  to 
assist  with  the  broad  array  of  semiconduetor  technologies. 

B,  ATLAS®/DECKBUILD 

The  ability  to  aceurately  simulate  a  semieonductor  deviee  is  eritical  to  industry 
and  research  environments.  The  ATLAS  deviee  simulator  is  specifically  designed  for 
2D  and  3D  modeling  to  inelude  eleetrical,  optieal  and  thermal  properties  within  a 
semiconductor  device.  ATLAS  provides  an  integrated  physies-based  platform  to 
analyze  DC,  AC  and  time-domain  responses  for  all  semieonductor-based  technologies. 
The  powerful  input  syntax  allows  the  user  to  design  any  semiconductor  device  using  both 
standard  and  user-defined  material  of  any  size  and  dimension.  ATLAS  also  offers  a 
number  of  useful  deviee  examples  to  assist  in  your  own  unique  design. 

DeckBuild  is  the  run  time  environment  used  to  input  a  command  file  or  deck  and 

(G) 

is  given  the  extension  “.in”.  To  run  ATLAS  in  the  DeckBuild  environment,  the  user 
must  first  call  the  ATLAS  simulator  with  the  command; 


go  atlas 


There  is  other  simulation  software  available  that  ean  be  used  in  conjunetion  with 
DeckBuild,  such  as  ATHENA  and  DevEdit  (Ligure  39),  but  ATLAS®  is  the  most 
appropriate  environment  for  this  research. 
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(Run  Time  Environment) 


Figure  39.  Silvaco  Inc.’s  simulation  environment,  from  [10]. 


C.  SYNTAX 

Once  ATLAS  is  called  there  is  a  syntax  structure  that  must  be  followed  in  order 
for  ATLAS  to  execute  the  command  file  successfully.  Figure  40  is  a  list  of  primary 
group  and  statement  structure  specifications.  Although  there  are  a  few  input  exceptions, 
the  input  file  statements  follow  the  general  format  of: 


<STATEMENT>  <PARAMETER>=<VALUE> 
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Statements 


Group 


1.  structure  Specification 


MESH 

REGION 

ELECTRODE 

DOPING 


2.  Material  Models  Specification 


MATERIAL 

MODELS 

CONTACT 

INTERFACE 


3.  Numerical  Method  Selection  METHOD 


4.  Solution  Specification 


LOG 

SOLVE 

LOAD 

SAVE 


5.  Results  Analysis  EXTRACT 

TONYPLOT 


Figure  40.  ATLAS  command  groups  with  the  primary  statements  in  each  group, 

from  [11]. 

D.  STRUCTURE  SPECIFICATION 
1,  Mesh 

The  mesh  statement  is  used  to  define  the  structure  in  an  inverted  2D  or  3D 
Cartesian  grid.  The  x-axis  is  positive  from  left  to  right,  the  y-axis  is  negative  from  bottom 
to  top  and  the  z-axis  is  positive  from  the  plane  of  the  x-y  axis  into  the  workspace.  The 
reason  for  the  inverted  y-axis  is  that  the  manufacturing  coordinates  are  usually  described 
as  depth  below  the  surface.  All  coordinates  are  entered  in  microns  and  the  spacing  is 
used  to  refine  the  sharpness  and  accuracy  at  an  assigned  location.  ATLAS  produces  a 
series  of  triangles  to  form  the  mesh  based  on  the  user  input  parameters.  Figure  41  shows 
the  mesh  of  a  tunnel  junction  along  with  the  parameters  entered  in  ATLAS  . 
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Figure  4 1 .  Tunnel  junetion  mesh. 


This  particular  tunnel  junction  is  only  Ipm  thick  and  lOpm  wide.  The  depth  (z- 
direction)  has  a  default  value  of  1pm  but  was  defined  as  1  x  10^  to  produce  an  output 
equal  to  one  square  centimeter,  which  corresponds  to  published  data.  Building  this 
semiconductor  as  a  1cm  x  1cm  square  would  only  add  nodes  to  the  structure  and 
ultimately  increase  the  calculation  time.  The  implementation  of  the  junction  to  a  solar 
cell  required  a  smaller  node  set  which  I  achieved  by  limiting  the  x-y  dimensions  to 
lOcm^ . 

The  accuracy  of  any  calculated  value  in  ATLAS  is  dependent  on  the  refinement 
of  the  mesh  at  the  junction  of  interest.  The  x.mesh  is  much  larger  than  the  y.mesh 
because  the  variables  being  measured/calculated  are  across  the  junction,  which  was 
defined  horizontally. 
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2,  Region 

The  region  statement  is  used  to  separate  the  initial  mesh  statement  into  distinet 
blocks  and  sets  the  initial  material  parameters  that  can  be  referred  to  later  by  region 
number.  All  meshed  areas  of  a  structure  must  be  assigned  a  region  and  the  regions  must 
be  ordered  from  lowest  to  highest  region.  For  instance,  region  5  cannot  be  defined  before 
region  4.  The  tunnel  junction  is  described  by  the  DeckBuild  commands  along  with  the 
region  lines  in  Figure  42. 


Data  from  tjsi.str 


Microns 

Figure  42.  Tunnel  junction  region. 

3,  Electrode 

Once  the  regions  are  set,  the  electrodes  must  be  assigned  to  the  desired  region  so 
that  it  can  be  electrically  analyzed.  The  electrode  can  be  assigned  to  any  region  or 
portion  of  a  region  including  the  “top”  or  “bottom”  as  it  is  in  this  research.  Because  this 
tunnel  junction  is  intended  to  be  used  as  an  interconnect  between  two  solar  cells,  it  was 
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unnecessary  to  assign  the  electrodes  to  metal  contacts  or  large  areas  that  would  constitute 
superfluous  calculations  or  error.  Therefore,  the  anode  was  assigned  to  the  top  of  the 
junction  and  cathode  was  assigned  to  the  bottom.  The  following  statements  were  used  to 
define  these  parameters: 


#electrodes 

electrode  name=anode  top 
electrode  name=cathode  bottom 


4.  Doping 

The  last  required  input  of  the  structure  specification  is  the  doping  statement.  The 
doping  statement  is  used  to  assign  the  doping  level  within  the  previously  assigned 
regions.  Various  properties  can  be  appended  to  the  doping  statement  to  specify  how  the 
semiconductor  was  doped  and  of  whether  the  region  is  n-  or  p-type.  The  desired  I-V 
relationship  in  forward  bias  will  determine  the  doping  characteristics  of  the 
semiconductor.  However,  the  following  statements  were  developed  for  a  Si  junction  with 
a  forward  bias  current  of  approximately  12A  at  0.06V: 

doping  uniform  p.type  conc=5el9  y.max=0.5 
doping  uniform  n.type  conc=le20  y.min=0.5 

Here  the  regions  were  already  defined  in  a  previous  statement  and  were  therefore 
referred  to  by  their  region  number.  Region  1  was  uniformly  doped  with  a  p-type  material 
with  a  concentration  of  5el9.  Similarly,  region  2  was  uniformly  doped  with  an  n-type 
material  at  a  concentration  of  le20. 

With  the  doping  statement  complete,  the  structure  specifications  are  complete  and 
it  is  possible  to  continue  with  the  simulation  following  the  guidelines  in  Figure  40. 
Unlike  the  structure  specifications,  the  remaining  four  groups  to  be  specified  do  not 
require  all  statements  to  run  without  error. 
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E.  QUANTUM  TUNNELING  MESH 


In  order  to  calculate  the  quantum  tunneling  through  a  potential  barrier,  finer  slices 
of  the  junction  must  be  taken.  ATLAS  does  this  with  the  definition  of  the  qtx.mesh  and 
qty.mesh  statements.  The  direction  of  tunneling  must  also  be  specified  with  a  different 
parameter  in  the  models  statement.  Figure  43  illustrates  the  extremely  refined  qty.mesh 
along  with  the  input  parameters  to  achieve  this  level  of  precision.  Although  the  qtx.mesh 
was  defined,  the  slices  are  only  extremely  refined  in  the  y-direction  because  this  is  the 
direction  of  tunneling.  The  new  rectangular  mesh  is  superimposed  over  the  original 
triangular  mesh  and  both  mesh  statements  will  be  evaluated  and  accounted  for  in  the 
tunneling  current.  As  with  all  ATLAS  spacing  parameters,  the  qty.mesh  spacing  is 
specified  as  pm. 
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qtx.mesh  loc=0  spac=l 

qtx. mesh  ioc=10  spac=l 

qty. mesh  iocs0.3S  spacsO.Oi 

qty.mesh  ioc=0.5  spac=0.0005 

qty.mesh  ioc=0.65  spac=0.0i 


Figure  43.  Tunnel  junction  qtx.mesh  and  qty.mesh. 
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F.  MATERIAL  MODELS  SPECIFICATION 

1,  Material 

The  material  statement  relates  physical  parameters  with  the  materials  assigned  to 
the  mesh.  The  important  material  parameters  for  most  standard  semiconductors  are 
already  defined  by  ATLAS  and  therefore  do  not  require  any  changes.  We  utilized  the 
extensive  model  libraries  in  ATLAS  since  the  developed  tunnel  junctions  were  based  in 
Si,  Ge,  and  GaAs  which  are  all  very  well  researched  and  documented. 

2,  Models 

The  models  statement  is  essential  to  the  accurate  modeling  of  a  particular 
phenomenon  because  it  sets  flags  for  ATLAS  to  indicate  the  inclusion  of  different 
mathematical  models,  physical  mechanisms  and  other  global  parameters  such  as  substrate 
temperature.  The  model  statement  and  parameters  used  in  this  research  are  listed  below: 

models  temperature=300  srh  fermi  ni . fermi  print  qtunn.dir=0  \ 
bbt. nonlocal  bbt. forward  bbt.nlderivs 


(C) 

As  with  all  ATLAS  statements,  a  complete  list  of  modeling  parameters  and  their 
definitions  can  be  found  in  the  ATLAS  User’s  Manual  [1 1].  The  most  notable 
parameters  for  the  calculation  of  the  quantum  tunneling  effect  are  qtunn.dir,  bbt.nonlocal, 
bbt.forward  and  bbt.nlderivs.  qtunn.dir  specifies  the  direction  of  tunneling  across  the 
junction.  Since  this  device  was  set  up  for  tunneling  across  the  y-axis,  the  parameter  value 
was  set  to  zero  (to  specify  tunneling  in  x-direction,  qtunn.dir  =  1).  The  bbt.nonlocal, 
bbt.forward  and  bbt.nlderivs  parameters  enable  the  non-local  band-to-band  tunneling 
model  to  be  used  in  forward  bias  and  allows  the  use  of  the  non-local  derivatives  in  the 
Jacobian  matrix  which  helps  eliminate  convergence  problems. 

3.  Contact 

Since  the  tunnel  junction  is  meant  to  be  used  as  an  interconnect  of  a  single  crystal 
multi-junction  solar  cell  there  was  no  need  to  specify  any  contact  parameters.  The 
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electrode  statements  assigning  the  cathode  and  anode  to  the  semiconductor  were 
sufficient  for  calculating  the  current  through  the  junction. 


G.  NUMERICAL  METHOD  SELECTION 

The  method  statement  specifies  the  numerical  methods  to  be  used  in  the 
calculations  specified  in  the  model  statement.  The  quantum  tunneling  through  a  junction 
is  highly  dependent  on  the  carrier  concentration.  Therefore,  it  is  very  important  to 
accurately  measure  the  carrier  concentration  for  each  iteration.  This  is  done  by  defining 
the  climit  parameter  in  the  method  statement,  which  specifies  the  normalization  factor  to 
limit  the  step  size  of  the  potential  update.  The  parameter  dvmax  sets  the  maximum 
allowed  potential  updates  per  iteration.  Reducing  dvmax  dampens  the  oscillations  in  the 
small  I-V  calculations. 

H.  SOLUTION  SPECIFICATION 

The  log/solve/save  statements  are  used  to  create  data  files  in  ATLAS 
simulations.  These  statements  work  together  to  be  provide  data  to  be  analyzed  by  other 
functions.  This  research  I  have  used  these  three  statements  as  a  tool  to  check  bias  points, 
energy  band  diagrams  and  carrier  concentrations  as  a  result  of  varying  stimuli. 

1.  Log 

The  log  statement  allows  all  terminal  characteristics  generated  by  a  solve 
statement  to  be  saved  to  a  file. 

2,  Solve 

The  solve  statement  specifies  which  bias  points  are  to  be  applied  to  produce  an 
output.  The  bias  points  can  be  set  in  a  number  of  different  way  including  step,  initial  and 
lambda  depending  on  what  stimulus  is  desired. 
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3. 


Save 


The  save  statement  is  used  to  save  all  node  point  data  into  an  output  file.  A 
common  use  of  the  log/solve/save  statements  for  the  tunnel  junction  studied  is  shown 
below: 


log  outfile=t j si . log 

solve  name=anode  vanode  =-l  vstep=0.1  vfinal=0.0 
solve  name=anode  vanode  =0.0  vstep=0.01  vfinal=0.2 
solve  name=anode  vanode  =0.2  vstep=0.1  vfinal=1.5 
save  outf ile=t j si . log 


Here  the  output  file  tjsi.log  has  been  set  aside  to  store  data  from  the  solve 
statements  and  then  saved  so  that  the  data  can  be  called  at  a  later  time.  The  solve 
statements  simply  sweep  the  voltage  from  -IV  to  1.5V  in  varying  voltage  steps 
depending  on  the  degree  of  accuracy  at  the  desired  range.  With  the  date  stored  in  an 
outlile,  it  is  ready  to  be  displayed  so  that  it  can  be  analyzed. 

I.  RESULTS  ANALYSIS 

The  tonyplot  statement  is  used  to  start  the  graphical  post-processor  tool.  Calling 
tonyplot  to  visualize  tjsi.log  displays  the  output  fde  in  Figure  44.  The  complete  tunnel 
junction  input  deck  can  be  found  in  Appendix  A. 
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Figure  44.  Silicon  tunnel  junction  I-V  curve  in  ATLAS  . 
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VI.  TUNNEL  JUNCTION  RESULTS 


After  the  successful  modeling  of  the  basic  forward  bias  tunneling  current  in  Si,  it 
was  necessary  to  test  the  validity  of  the  model.  This  was  accomplished  by  measuring  the 
response  of  different  materials  and  doping  levels  to  similar  bias  conditions. 

A.  SEMICONDUCTOR  MATERIAL 

As  an  interconnect  between  two  dissimilar  solar  cells  it  is  important  that  a  tunnel 
junction  can  be  modeled  for  materials  other  than  Si.  Typical  tunnel  junctions  used  today 
are  GaAs  and  InGaP  due  to  their  close  lattice  matching.  The  ATLAS  input  deck 
presented  here  was  developed  so  that  cell  parameters  are  easily  changed  and  the  tonyplot 
outputs  are  easy  to  compare. 

The  doping  levels’  initial  calculations  were  based  on  the  Fermi-Dirac  equations 
presented  in  Section  II.  Figures  45  and  46  illustrate  the  tunnel  junction  energy  band 
diagrams  along  with  their  associated  I-V  curves  for  GaAs  and  InGaP  respectively. 


0.1  0.2  0.3  0.4  0.6  0.6  0.7  0.8  0.9  1 


0  0.2  0.4  0.6  0.8  1  1.2 


Distance  along  line 


Anode  Voltage  (V) 


3  3 

Figure  45.  GaAs  tunnel  junction  (p  =  3 e20/cm  n  =  9el8/cm) 

{PJcvcv  versus  volts  in  a  Ipm  thick  junction). 
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Figure  46.  InGaP  tunnel  junction  (p  =  2e20/cm^  n  =  9el8/cm^) 
{PJcvcv  versus  volts  in  a  Ipm  thick  junction). 


By  fine-tuning  the  doping  levels,  an  output  current  of  approximately  200mA  was 
realized  for  both  junctions.  Though  this  current  is  slightly  high  for  a  typical  current 
produced  by  a  solar  cell,  we  simply  solved  for  200mA  as  a  baseline.  With  the  successful 
modeling  of  the  two  junctions,  we  were  able  to  manipulate  the  doping  levels  and  analyze 
the  different  outputs. 

B.  TUNNEL  JUNCTION  DOPING 

The  doping  levels  of  the  tunnel  junctions  were  changed  for  both  the  n-  and  p- 
regions  to  shift  the  Fermi  level  and  ultimately  vary  the  response  to  the  given  bias 
condition.  Figures  47  through  52  illustrate  the  reaction  to  increased  doping  of  the  GaAs 
junction.  Figures  53  through  58  show  similar  results  for  an  InGaP  junction. 

Initially,  the  p-region  was  varied  from  a  doping  level  within  3kT  (~0.077eV)  of 
the  valence  band,  then  increased  by  a  reasonable  factor  (between  2  and  3)  to  show  the 
response  in  the  I-V  curve.  Pushing  the  Fermi  level  further  into  the  valence  band 
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corresponds  to  an  increased  number  of  available  holes  that  can  recombine  with  electrons 
tunneling  through  the  junction  and  therefore  increases  the  current. 

The  n-region  was  varied  in  a  similar  manner.  However,  pushing  the  Fermi  level 
further  into  the  conduction  band  increases  the  number  of  electrons  available  to  tunnel 
through  the  junction  and  similarly  increases  the  current. 


1,  GaAs  Tunnel  Junction 


a,  p-region  variation: 


Figure  47. 


2 

GaAs  IV  characteristics  and  energy  band  diagram  (p  =  le20/cm  n  = 
9el8/cm^)  {AJcvcv  versus  volts  in  a  1pm  thick  junction). 
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Figure  48. 
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GaAs  IV  characteristics  and  energy  band  diagram  (p  =  3e20/cm  n 
9el8/cm^)  {AJcvcv  versus  volts  in  a  1pm  thick  junction). 
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Figure  49. 


GaAs  IV  characteristics  and  energy  band  diagram  (p  =  6e20/cm  n 
9el8/cm^)  (A/cm^  versus  volts  in  a  1pm  thick  junction). 


68 


b, 


n-region  variation 
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Figure  50. 


GaAs  IV  characteristics  and  energy  band  diagram  (p  =  3e20/cm^  n 
6el8/cm^)  {AJcvcv  versus  volts  in  a  1pm  thick  junction). 


Figure  5 1 .  GaAs  IV  characteristics  and  energy  band  diagram  (p  =  3e20/cm  n 

9el8/cm^)  {AJcvcv  versus  volts  in  a  1pm  thick  junction). 
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Figure  52. 


2. 


Figure  53. 
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GaAs  IV  charaeteristies  and  energy  band  diagram  (p=3e20/cm  n=2el9/em  ) 
(A/em2  versus  volts  in  a  Igm  thiek  junetion). 


InGaP  Tunnel  Junction 
a,  p-region  variation: 
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InGaP  IV  eharacteristies  and  energy  band  diagram  (p  =  8el9/em  n  = 

3  2 

9el8/em  )  (A/cm  versus  volts  in  a  Ipm  thick  junction). 
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Figure  54.  InGaP  IV  characteristics  and  energy  band  diagram  (p  =  2e20/cm  n 

3  2 

9el8/cm  )  (A/cm  versus  volts  in  a  Ipm  thick  junction). 


Figure  55. 


InGaP  IV  characteristics  and  energy  band  diagram  (p  =  4e20/cm  n 
9el8/cm^)  {AJcm  versus  volts  in  a  Ipm  thick  junction). 
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b, 


n-region  variation 


Figure  56.  InGaP  IV  characteristics  and  energy  band  diagram  (p  =  2e20/cm  n 
6el8/cm^)  {AJcvcv  versus  volts  in  a  1pm  thick  junction). 
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Figure  57. 


InGaP  IV  characteristics  and  energy  band  diagram  (p  =  2e20/cm  n 

3  2 

9el8/cm  )  (A/cm  versus  volts  in  a  1pm  thick  junction). 
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Figure  58.  InGaP  IV  characteristics  and  energy  band  diagram  (p  =  2e20/cm  n  = 
2el9/cm^)  {AJcvcv  versus  volts  in  a  1pm  thick  junction). 


C.  SUMMARY 

With  the  characteristic  tunnel  junction  curve  successfully  modeled,  I  chose  to  take 
my  research  one  step  further  and  implement  the  tunnel  junction  into  a  solar  cell 
simulation  to  verify  the  operation  which  will  be  discussed  in  the  next  chapter. 
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VII.  SOLAR  CELL  IMPLEMENTATION 


The  solar  cell  implementation  was  based  on  research  completed  by  Panayiotis 
Michalopoulos  [2]  in  2002.  Michalopoulos  was  a  student  of  Prof.  Sherif  Michael  at  the 
Naval  Postgraduate  School  and  his  research  was  developed  using  Silvaco  software. 

A,  SOLAR  RADIATION 

The  sun’s  energy  is  derived  from  nuclear  fusion  reactions  where  6  x  lO'^kg  of 
hydrogen  (H2)  is  converted  to  helium  (He)  every  second.  The  mass  lost  due  to  the 
conversion  is  approximately  4  x  10  kg.  The  energy  produced  by  this  loss  can  be 
determined  by  Einstein’s  famous  mass-energy  equation  (E=mc  )  and  is  found  to  be  4  x 
10  J.  This  energy  is  emitted  as  electromagnetic  radiation  primarily  at  wavelengths  from 
the  ultraviolet  to  infrared  and  radio  spectra  (0.2  -  3pm)  [12]. 

The  intensity  of  the  sun’s  energy  depends  on  the  location  of  where  it  is  measured. 

2 

Outside  the  earth’s  atmosphere,  the  intensity  is  measured  to  be  approximately  1365W/m 
which  produces  a  specific  spectral  response  referred  to  as  air  mass  zero  (AMO).  As  the 
energy  transits  the  earth’s  atmosphere,  it  is  attenuated  by  ultraviolet  absorption  due  to 
ozone,  scattering  due  to  airborne  dust,  and  infrared  absorption  due  to  water  vapor.  AMI 
and  AM  1.5  is  used  for  terrestrial  applications  with  incident  solar  power  ranging  from  925 
to  963W/ml 

B,  I-V  CURVE 

As  discussed  in  Chapter  II,  an  electron  must  release  energy  to  go  from  an  outer 
shell  (higher  energy  state)  to  an  inner  shell  (lower  energy  state).  Conversely,  if  a  certain 
level  of  energy  is  absorbed,  such  as  solar  energy,  an  electron  may  transit  from  an  inner 
shell  to  an  outer  shell  or  valence  shell.  In  an  intrinsic  semiconductor  this  creation  of 
carriers  leads  to  an  increase  in  conductivity.  However,  when  a  p-n  junction  is  formed  in  a 
semiconductor  the  carriers  will  separate  to  create  electron  hole  pairs  (EHPs).  The 
minority  carrier  will  be  swept  across  the  junction  due  to  the  electrostatic  field  formed  by 
the  depletion  region  and  the  excess  majority  carriers  created  by  the  photogeneration  will 
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develop  a  voltage  differential  between  the  two  regions.  If  the  junction  is  then  connected 
to  a  load,  current  will  flow  (photocurrent).  This  phenomenon  is  known  as  photovoltaics 
because  of  the  voltage  potential  generated  by  the  photon  (solar)  energy. 

As  photons  generate  EHPs  in  a  p-n  junction,  the  standard  diode  dark  current 
characteristic  (Figure  59a-blue)  is  shifted  down  to  produce  the  illuminated  current 
characteristic  (Figure  59a-red).  Because  a  solar  cell  is  usually  presented  as  a  power 
supply,  the  illuminated  characteristic  is  most  often  reflected  across  the  x-axis  to  show  that 
the  photocurrent  conventionally  flows  in  the  opposite  direction  of  the  standard  diode 
current  (Figure  7.1b). 


Figure  59.  (a)  Dark  (blue)  and  illuminated  (red)  current  characteristic  curve  of  the  solar 

cell,  (b)  Solar  cell  current  characteristic  curve  with  FF. 


The  number  of  FHPs  generated  in  the  semiconductor  is  proportional  to  the 
photons  entering  the  semiconductor.  However,  many  of  the  FHPs  recombine  inside  the 
material  before  it  is  collected  and  are  therefore  not  used  to  produce  photocurrent.  Thus, 
the  maximum  obtained  power  (Pobtamed)  than  the  maximum  theoretical  power 
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(^theoiy  )•  The  ratio  of  obtained  power  to  theoretieal  power  is  ealled  fill  faetor  or  FF  (Eqn. 
7.1).  The  FF  ean  be  seen  graphieally  in  Figure  59  (b). 


FF  = 


^  obtained 


^  theory 


I  V 

max  max 

I  V 

^sc  oc 


Eqn.  7.1 


The  power  eonversion  effieieney  (n)  of  a  solar  eell  is  given  as  the  ratio  of  the 
maximum  realized  power  to  the  ineident  power  ( ). 


P  I  V  EE 

n  =  max  _  ^SC  ^ 

P„„  P„„ 


Eqn.  7.2 


To  maximize  the  power  eonversion  effieieney,  the  EE,  Vqc  and  fie  must  be 
optimized. 

C.  SINGLE  CELL 

Initially  a  single  eell  for  eaeh  material  in  the  multi-junetion  eell  was  modeled  in 
ATLAS  .  This  simulation  was  based  on  the  eode  already  developed  for  a  Naval 
Postgraduate  Sehool  thesis  by  Panayiotis  Miehalopoulos  in  Mareh  2002.  Both  eells  were 
illuminated  at  AMO.  The  doping  levels  and  thieknesses  are  illustrated  in  Eigure  60  and 
62  followed  by  their  eorresponding  I-V  eurve  (Eigures  61  and  63). 
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Figure  60.  Optimized  InGaP  cell  structure. 
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Figure  61.  Optimized  InGaP  I-V  curve.  (A/cm  versus  voltage). 
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Figure  62.  Optimized  GaAs  eell  structure. 
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Figure  63.  Optimized  GaAs  I-V  curve.  (A/cm  versus  voltage). 
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The  InGaP  cell  achieved  a  Voc  =  1.43V  and  an  Isc  =  19niA/cm^.  The  GaAs  cell 

2 

achieved  a  Voc  =  I.IV  and  an  Isc  =  27.4niA/cm  .  Both  cells  were  successfully  modeled 
with  5%  higher  power  than  published  data  and  2%  higher  power  than  the  results  achieved 
by  Michalopoulos  [2], 

With  the  single  cells  successfully  modeled,  the  tunnel  junction  was  implemented 
into  the  ATLAS®  code. 

D.  MULTI  CELLS 

Similar  to  the  single  solar  cells,  the  multi-junction  cell  in  Figure  64  was 
illuminated  at  AMO  and  the  I-V  curve  was  recorded  (Figure  65). 
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Figure  64.  Multi-junction  solar  cell  (InGaP/GaAs)  with  tunnel  junction. 
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Figure  65.  Multi-junction  solar  cell  (InGaP/GaAs)  I-V  curve.  (A/cm  versus  voltage). 


The  multi-junction  solar  cell  model  achieved  a  Voc  =  2.5V  and  an  Isc  = 
18.7mA/cm  .  As  expected,  the  voltage  is  the  sum  of  the  two  cells  in  the  single  cell 
models  accounting  for  the  small  voltage  drop  across  the  tunnel  junction.  Ideally,  both 
cells  would  generate  the  same  current  as  in  the  single  cell  model  and  the  realized  current 
in  the  multi-junction  cell  would  be  the  lower  of  the  two.  However,  because  the  higher 
energy  photons  are  absorbed  by  the  InGaP  top-cell,  the  GaAs  bottom-cell  is  not  able  to 
generate  as  many  EHPs  and  therefore  realizes  a  slightly  smaller  current  than  ideal. 

Taking  the  research  one  step  further,  the  tunnel  junction  was  implemented  in  the 
optimized  multi-junction  cell  model.  This  was  done  so  that  the  modeled  output  could  be 
compared  to  empirical  data  in  [3].  The  optimized  cell  has  an  antireflective  coating  to 
absorb  a  greater  number  of  photons  and  a  back  surface  reflector  to  reflect  some  of  the 
unutilized  photons  back  into  the  semiconductor  material.  These  additions  help  to  create 
more  EHPs  and  ultimately  a  higher  current  and  power  density. 

Experimental  data  for  an  InGaP/GaAs  photovoltaic  cell  produced  a  Voc  =  2.488V 
and  an  Isc  =  23mA/cm  [3].  The  optimized  multi-junction  cell  model  achieved  a  Voc  = 
2.5V  and  an  Isc  =  22mA/cm^  with  a  EE  =  84%  and  rj  =  29%  (Eigure  66).  This  is  only 
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3.8%  greater  power  density  than  the  experimental  data  and  demonstrates  the  ability  to 
model  a  tunnel  junetion  for  use  in  photovoltaie  optimization. 
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Figure  66. 


Data  from  multilV.log 


Anode  Int.  Voltage  (V) 


Optimized  multi-junetion  solar  eell  (InGaP/GaAs)  I-V  eurve. 
(A/em  versus  voltage) 
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VIII.  RECOMMENDATIONS  AND  CONCLUSION 


A.  RECOMMENDATIONS 

1,  Physical  Experiment 

The  tunnel  junetion  model  presented  here  has  produeed  very  promising  results. 
For  this  research  the  most  important  characteristic  of  the  tunnel  junction  is  the  forward 
bias  electron  tunneling  which  was  successfully  modeled  using  Silvaco  Inc’s  ATLAS 
software.  Unfortunately,  time  and  resources  did  not  allow  for  a  physical  model 
comparison.  Although  an  implemented  junction  was  compared  to  empirical  values,  the 
tunnel  itself  was  only  compared  to  theoretical  values  and  responses  based  on 
mathematical  models. 

After  further  refinement  and  modeling,  it  is  recommended  that  a  project  be  funded 
to  produce  a  junction  with  the  precise  specifications  so  that  it  may  be  tested  and 
compared  to  the  modeled  parameters.  Developing  a  more  efficient  tunnel  junction  and 
continuing  to  test  the  accuracy  of  the  models,  will  lead  to  the  development  of  more 
efficient  solar  cells. 

2.  Other  Software 

Throughout  this  study,  a  number  of  competing  software  programs  were  used  in  an 
effort  to  model  the  forward  bias  response  of  a  tunnel  junction.  Some  were  geared  more 
toward  other  aspects  of  the  semiconductor  industries  such  as  MEMs  design  or  computer 
architecture  and  were  therefore  difficult  to  use  to  implement  a  tunnel  junction  into  a  solar 
cell.  Other  programs  were  solar  cell  centric  yet  did  not  have  the  capacity  or  tools  to 
realize  the  quantum  tunneling  effect  of  the  degenerately  doped  junction. 

There  are  other  software  packages  available  that  are  comparable  to  Silvaco’s 
ATLAS  and  would  be  beneficial  to  this  research.  It  is  recommended  that  the  utility  of 
these  software  packages  be  explored  and  compared. 
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3.  Temperature  and  Radiation  Effects 

The  use  of  photovoltaies  in  spaee  applieations  is  very  important.  Therefore,  the 
tunnel  junetion  should  be  adapted  and  tested  under  extreme  eonditions  to  ensure  the 
aeeuraey  of  the  model.  The  effeets  of  temperature  and  radiation  on  a  tunnel  junetion 
should  be  explored  to  ensure  no  anomalies  exist. 

B,  CONCLUSION 

The  value  of  using  simulation  software  sueh  as  ATLAS®  to  model 
semieonduetors  has  revolutionized  the  semieonduetor  industry.  Not  only  are  time  and 
money  saved  by  minimizing  the  need  to  manufaeture  and  test  ehanges  to  a  design,  but 
effleieney  is  being  pushed  to  the  limit. 

It  has  been  have  shown  throughout  this  thesis  that  it  is  possible  to  model  sueh 
effeets  as  quantum  tunneling  in  a  photovoltaie  solar  eell.  To  this  point,  the  tunnel 
junetion  in  a  multi-junetion  eell  was  modeled  as  a  vaeuum  to  insulate  the  layers,  then 
hardwired  to  pass  the  eurrent  between  eells.  Although  this  is  effeetive,  the  eell  being 
modeled  is  not  fully  optimized. 

A  paper  entitled  Multi- Junction  Photovoltaic  Model  Optimizations  for  Space  and 
Solar  Concentrator  Applications”  [13],  was  presented  by  Professor  Miehael  and  the 
author  on  September  3,  2008,  at  the  23rd  European  Photovoltaie  Solar  Energy 
Conferenee,  that  was  held  in  Valeneia,  Spain.  The  paper  foeused  on  the  applieation  of 
multi-junetion  photovoltaie  eells  using  tunnel  junetion  teehnology  for  use  in  spaee 
applieations.  The  tunnel  junetion  modeling  has  the  potential  to  push  photovoltaie 
teehnology  to  a  new  level  of  optimization  and  power  density  and  it  is  neeessary  to  model 
this  phenomenon  aeeurately. 
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APPENDIX.  ATLAS®  SOURCE  CODE 


A.  MAIN  STRUCTURE 

go  atlas 

#  Definition  of  constants 

#  Mesh 

#  X-Mesh 

#  Y-Mesh 

#  Regions 

#  Electrodes 

#  Doping 

#  QT  Mesh 

#  Material  properties 

#  Models 

#  Light  beams 

#  Solving 


B,  TUNNEL  JUNCTION 

go  atlas 

1.  X-Mesh 

mesh  width=le7 
#x-mesh 

x.mesh  loc=0  spac=l 

x. mesh  loc=10  spac=l 

2,  Y-Mesh 

y. mesh  loc=0.0  spac=0.01 

y.mesh  loc=0.25  spac=0.001 

y.mesh  loc=0.75  spac=0.001 

y.mesh  loc=l . 0  spac=0.01 

3.  Regions 

region  num=l  material=InGaP  x.min=0  x.max=10  y.min=0.5  y.max=l 
region  num=2  material=InGaP  x.min=0  x.max=10  y.min=0  y.max=0.5 

4,  Electrodes 

elec  name=anode  top 

elec  name=cathode  bottom 

doping  uniform  p.type  conc=2e20  y.max=0.5 
doping  uniform  n.type  conc=2el9  y.min=0.5 
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5. 


QT-Mesh 


qtx.mesh  loc=0  spac=l 
qtx.mesh  loc=10  spac=l 


qty.mesh  loc=0.35 
qty.mesh  loc=0.5 
qty.mesh  loc=0.65 


spac=0 . 01 
spac=0 . 0005 
spac=0 . 01 


6.  Models/Output 


models  temperature=300  srh  fermi  ni.fermi  print  \ 
qtunn.dir=0  bbt. nonlocal  bbt. forward  bbt.nlderivs 

output  band. temp  traps  u.srh  taurn  taurp  band.param  con. band  val.band  \ 
qfn  qfp 

log  outf ile=t j InGaP . str 
solve  init 

save  outf ile=tj InGaP . str 
tonyplot  tj InGaP. str 


method  climit=1.0  dvmax=0.05 


solve  init 
solve  prev 

log  outf ile=tj InGaP . log 

7.  Solve 

solve  name=anode  vanode 
solve  name=anode  vanode 
solve  name=anode  vanode 
tonyplot  tjInGaP.log 


=-l  vstep=0.1  vfinal=0.0 
=0.0  vstep=0.01  vfinal=0.2 
=0.2  vstep=0.1  vfinal=l.l 


8.  Plot 

save  outf ile=tj InGaP . str 
tonyplot  tj InGaP. str 


C.  COMMON  SECTIONS 
1.  X-Mesh 

mesh  space. mult=l 

#  X-Mesh:  surface=500  um2  =  1/200,000  cm2 
x.mesh  loc=-250  spac=50 
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x.mesh  loc=0  spac=10 
x.mesh  loc=250  spac=50 


2,  Material  Properties 

material  TAUN=le-7  TAUP=le-7  COPT=1.5e-10  AUGN=8.3e-32  AUGP=1.8e-31 

#  Vacuum 

material  material=Vacuum  real . index=3 . 3  imag . index=0 

#  Ge 

material  material=Ge  EG300=0.67  PERMITTIVITY=1 6  AFFINITY=4 
material  material=Ge  MUN=3900  MUP=1800 
material  material=Ge  NC300=1 . 04el9  NV300=6el8 
material  material=Ge  index . file=Ge . opt 

#  GaAs 

material  material=GaAs  EG300=1.42  PERMITTIVITY=13 . 1  AFFINITY=4 . 07 
material  material=GaAs  MUN=8800  MUP=400 
material  material=GaAs  NC300=4.7el7  NV300=7el8 
material  material=GaAs  index . file=GaAs . opt 

#  InGaP 

material  material=InGaP  EG300=1.9  PERMITTIVITY=1 1 . 62  AFFINITY=4 . 1 6 
material  material=InGaP  MUN=1945  MUP=141 
material  material=InGaP  NC300=1.3e20  NV300=1 . 28el9 
material  material=InGaP  index . f ile=InGaP-l . 9 . opt 

#  AllnP  (=InAsP) 

material  material=InAsP  EG300=2.4  PERMITTIVITY=1 1 . 7  AFFINITY=4.2 
material  material=InAsP  MUN=2291  MUP=142 
material  material=InAsP  NC300=1 . 08e20  NV300=1 . 28el9 
material  material=InAsP  index . file=AlInP . opt 

#  AlInGaP  (=InAlAsP) 

material  material=InAlAsP  EG300=2.4  PERMITTIVITY=1 1 . 7  AFFINITY=4.2 
material  material=InAlAsP  MUN=2150  MUP=141 
material  material=InAlAsP  NC300=1.2e20  NV300=1 . 28el9 
material  material=InAlAsP  index . file=AlInP . opt 

3,  Models 

models  term!  ni.fermi  qtunn.dir=0  bbt. nonlocal  bbt. forward  bbt.nlderivs 
bbt . std 

4,  Light  Beams 

beam  num=l  x.origin=0  y.origin=-5  angle=90  \ 

power . file=AM0silv . spec  wavel . start=0 . 2 1  wavel.end=4  wavel.num=50 
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D,  INGAAS  /  GAAS  SINGLE  AND  DUAL  JUNCTION  CELL 


1,  Bottom  Cell  (GaAs) 
a,  Y-Mesh 


#  Vacuum 

y.mesh  loc=-0.15  spac=0.001 

#  Window  (0.05  um) 
y.mesh  loc=-0.1  spac=0.01 

#  Emitter  (0.1  um) 
y.mesh  loc=0  spac=0.01 

#  Base  (3  um) 

y.mesh  loc=1.5  spac=0.3 
y.mesh  loc=3  spac=0.01 

#  BSE  (0.1  um) 

y.mesh  loc=3.1  spac=0.01 

#  Buffer  (0.3  um) 
y.mesh  loc=3.4  spac=0.05 

#  Substrate  (300  um) 
y.mesh  loc=303.4  spac=50 


b.  Regions 

#  Window  AllnP  (=InAsP) 

region  num=l  material=InAsP  x.min=-250  x.max=250  y.min=-0.15  y.max=-0.1 

#  Emitter 

region  num=2  material=GaAs  x.min=-250  x.max=250  y.min=-0.1  y.max=0 

#  Base 

region  num=3  material=GaAs  x.min=-250  x.max=250  y.min=0  y.max=3 

#  BSE 

region  num=4  material=InGaP  x.min=-250  x.max=250  y.min=3  y.max=3.1 

#  Buffer 

region  num=5  material=GaAs  x.min=-250  x.max=250  y.min=3.1  y.max=3.4 

#  Substrate 

region  num=6  material=GaAs  x.min=-250  x.max=250  y.min=3.4  y.max=303.4 

c.  Electrodes 

electrode  name=cathode  x.min=-250  x.max=250  y.min=-0.15  y.max=-0.15 
electrode  name=anode  x.min=-250  x.max=250  y.min=303.4  y.max=303.4 

d.  Doping 


#  Window 

doping  uniform  region=l  n.type 

#  Emitter 

doping  uniform  region=2  n.type 

#  Base 

doping  uniform  region=3  p.type 

#  BSE 

doping  uniform  region=4  p.type 


conc=lel 9 
conc=2el 8 
conc=lel7 
conc=2el 8 
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#  Buffer 

doping  uniform  region=5  p.type  conc=7el8 

#  Substrate 

doping  uniform  region=6  p.type  conc=lel9 

e.  Solving 

method  newton  maxtraps=10  itlimit=100  climit=l  dvmax=0.05 
solve  initial 
solve  bl=0 . 9 

log  outf=multiIV. log 

solve  vanode=0.0  name=anode  vstep=0.025  vfinal=2.7 
extract  init  inf ile="multiIV. log" 

extract  name="Jsc"  y.val  from  curve (v . "anode" ,  1. "anode") 
X . val=0 . 0 

extract  name="Voc"  x.val  from  curve (v . "anode" ,  1. "anode") 
y . val=0 . 0 

2,  Top  Cell  (InGaP) 
a.  Y-Mesh 


#  Vacuum 

y.mesh  loc=-0.87  spac=0.003 

#  Window  (0.03  um) 

y.mesh  loc=-0.84  spac=0.003 

#  Emitter  (0.05  um) 
y.mesh  loc=-0.79  spac=0.003 

#  Base  (0.55  um) 
y.mesh  loc=-0.5  spac=0 . 1 
y.mesh  loc=-0.24  spac=0.003 

#  BSE  (0.03  um) 

y.mesh  loc=-0.21  spac=0.003 

#  Buffer  (0.03  um) 

y.mesh  loc=-0.18  spac=0.002 


b.  Regions 

#  Window  AllnP  (=InAsP) 

region  num=l  material=InAsP  x.min=-250  x.max=250  y.min=-0 
0.84 

#  Emitter 

region  num=2  material=InGaP  x.min=-250  x.max=250  y.min=-0 
0.79 

#  Base 

region  num=3  material=InGaP  x.min=-250  x.max=250  y.min=-0 
0.24 

#  BSE 

region  num=4  material=InGaP  x.min=-250  x.max=250  y.min=-0 

0.21 


where 

where 


87  y.max= 

84  y.max= 

79  y.max= 

24  y.max= 
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#  Buffer  AllnP  (=InAsP) 

region  num=5  material=InAsP  x.min=-250  x.max=250  y.inin=-0.21  y.max=- 
0.18 

c.  Electrodes 

electrode  name=cathode  x.min=-250  x.max=250  y.min=-0.87  y.max=-0.87 
electrode  name=anode  x.min=-250  x.max=250  y.min=-0.18  y.max=-0.18 

d.  Doping 

#  Window 

doping  uniform  region=l  n.type  conc=1.95el8 

#  Emitter 

doping  uniform  region=2  n.type  conc=2el8 

#  Base 

doping  uniform  region=3  p.type  conc=1.5el7 

#  BSE 

doping  uniform  region=4  p.type  conc=2el8 

#  Buffer 

doping  uniform  region=5  p.type  conc=0.95el8 

e.  Solving 

(1)  IscandVoc 

method  newton  maxtraps=10  itlimit=100  climit=l  dvmax=0.05 
solve  initial 
solve  bl=0.9 

log  outf=multiIV. log 

solve  vanode=0.0  name=anode  vstep=0.025  vfinal=2.7 
extract  init  inf ile="multiIV. log" 

extract  name="Jsc"  y.val  from  curve (v . "anode" ,  i. "anode")  where 
X . val=0 . 0 

extract  name="Voc"  x.val  from  curve (v . "anode" ,  i. "anode")  where 
y . val=0 . 0 


3.  Stacked  Cell 

a.  Y-Mesh 

#  Vacuum 

y.mesh  loc=-0.87  spac=0.003 

#  Window  (0.03  um) 

y.mesh  loc=-0.84  spac=0.003 

#  Emitter  (0.05  um) 
y.mesh  loc=-0.79  spac=0.003 

#  Base  (0.55  um) 
y.mesh  loc=-0.5  spac=0.1 
y.mesh  loc=-0.24  spac=0.003 
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#  BSF  (0.03  urn) 

y.mesh  loc=-0.21  spac=0.003 

#  Buffer  (0.03  um) 

y.mesh  loc=-0.18  spac=0.002 

#  Vacuum  (0.015  um) 

y.mesh  loc=-0.165  spac=0.002 

#  Vacuum  (0.015  um) 
y.mesh  loc=-0.15  spac=0.001 

#  Window  (0.05  um) 
y.mesh  loc=-0.1  spac=0.01 

#  Emitter  (0.1  um) 
y.mesh  loc=0  spac=0.01 

#  Base  (3  um) 

y.mesh  loc=1.5  spac=0.3 
y.mesh  loc=3  spac=0.01 

#  BSF  (0.1  um) 

y.mesh  loc=3.1  spac=0.01 

#  Buffer  (0.3  um) 
y.mesh  loc=3.4  spac=0.05 

#  Substrate  (300  um) 
y.mesh  loc=303.4  spac=50 


b.  Regions 


#  Window  AllnP  (=InAsP) 

region  num=l  material=InAsP  x.min=-250  x.max=250  y.min=-0.87  y.max=- 
0.84 

#  Emitter 

region  num=2  material=InGaP  x.min=-250  x.max=250  y.min=-0.84  y.max=- 
0.79 

#  Base 

region  num=3  material=InGaP  x.min=-250  x.max=250  y.min=-0.79  y.max=- 
0.24 

#  BSF 


region  num=4  material=InGaP  x.min=-250  x.max=250  y.min=-0.24  y.max=- 

0.21 

#  Buffer  AllnP  (=InAsP) 

region  num=5  material=InAsP  x.min=-250  x.max=250  y.min=-0.21  y.max=- 
0.18 

#  Vacuum 

region  num=6  material=Vacuum  x.min=-250  x.max=250  y.min=-0.18 
y .max=-0 . 165 

region  num=7  material=Vacuum  x.min=-250  x.max=250  y.min=-0.165 
y .max=-0 . 15 

#  Window  AllnP  (=InAsP) 

region  num=8  material=InAsP  x.min=-250  x.max=250  y.min=-0.15 
y .max=-0 . 1 

#  Emitter 

region  num=9  material=GaAs  x.min=-250  x.max=250  y.min=-0.1  y.max=0 

#  Base 

region  num=10  material=GaAs  x.min=-250  x.max=250  y.min=0  y.max=3 

#  BSF 


region  num=ll  material=InGaP  x.min=-250  x.max=250  y.min=3  y.max=3.1 
#  Buffer 

region  num=12  material=GaAs  x.min=-250  x.max=250  y.min=3.1  y.max=3.4 
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#  Substrate 

region  num=13  material=GaAs  x.min=-250  x.max=250  y.inin=3.4  y.max=303.4 

c.  Electrodes 

electrode  name=cathode  x.min=-250  x.max=250  y.min=-0.87  y.max=-0.87 
electrode  name=cc  x.min=-250  x.max=250  y.min=-0.18  y.max=-0.18 
electrode  name=ee  x.min=-250  x.max=250  y.inin=-0.15  y.inax=-0.15 
electrode  name=anode  x.min=-250  x.max=250  y.min=303.4  y.max=303.4 

d.  Doping 

#  Window 

doping  uniform  region=l  n.type  conc=1.95el8 

#  Emitter 

doping  uniform  region=2  n.type  conc=2el8 

#  Base 

doping  uniform  region=3  p.type  conc=1.5el7 

#  BSE 

doping  uniform  region=4  p.type  conc=2el8 

#  Buffer 

doping  uniform  region=5  p.type  conc=0.95el8 

#  Window 

doping  uniform  region=8  n.type  conc=lel9 

#  Emitter 

doping  uniform  region=9  n.type  conc=2el8 

#  Base 

doping  uniform  region=10  p.type  conc=lel7 

#  BSE 

doping  uniform  region=ll  p.type  conc=2el8 

#  Buffer 

doping  uniform  region=12  p.type  conc=7el8 

#  Substrate 

doping  uniform  region=13  p.type  conc=lel9 

e.  Solving 

method  newton  maxtraps=10  itlimit=100  climit=l  dvmax=0.05 
solve  initial 
solve  bl=0.9 

log  outf=multiIV. log 

solve  vanode=0.0  name=anode  vstep=0.025  vfinal=2.7 
extract  init  inf ile="multiIV. log" 

extract  name="Jsc"  y.val  from  curve (v . "anode" ,  i. "anode")  where 
X . val=0 . 0 

extract  name="Voc"  x.val  from  curve (v . "anode" ,  i. "anode")  where 
y. val=0 . 0 
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4. 


Tunnel  Junction 


a.  Y-Mesh 

y.mesh  loc=-0.18  spac=0.002 

#  Tunnel  emitter  (0.015  um) 
y.mesh  loc=-0.165  spac=0.002 

#  Tunnel  base  (0.015  um) 
y.mesh  loc=-0.15  spac=0.001 

b.  Regions 

region  num=l  material=InGaP  x.min=-250  x.max=250  y.min=-0.18  y.max= 
0.165 

region  num=2  material=InGaP  x.min=-250  x.max=250  y.min=-0.165  y.max 
0.15 


c.  Electrodes 


electrode  name=cathode  x.min=-250  x.max=250  y.min=-0.18  y.max=-0.18 
electrode  name=anode  x.min=-250  x.max=250  y.min=-0.15  y.max=-0.15 

d.  Doping 

doping  uniform  region=l  p.type  conc=8el8 
doping  uniform  region=2  n.type  conc=lel9 


e,  QT  Mesh 


qtx .mesh 
qtx .mesh 
qtx .mesh 


loc=-250 

loc=0 

loc=250 


spac=50 

spac=10 

spac=50 


qty.mesh  loc==-0.165 
qty.mesh  loc=-0.1575 
qty.mesh  loc=-0.15 


spac= .001 
spac= .0005 
spac= .001 


f.  Solving 

method  newton  maxtraps=10  itlimit=100  climit=l  dvmax=0.05 
solve  initial 
solve  bl=0.9 

log  outf=multiIV. log 

solve  vanode=0 . 0  name=anode  vstep=0.025  vfinal=2.7 
extract  init  inf ile="multiIV. log" 

extract  name="Jsc"  y.val  from  curve (v . "anode" ,  1. "anode")  where 
X . val=0 . 0 

extract  name="Voc"  x.val  from  curve (v . "anode" ,  i. "anode")  where 
y . val=0 . 0 
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5. 


MJ  Cell 


a.  Y-Mesh 


#  Vacuum 

y.mesh  loc=-0.87  spac=0.003 

#  Window  (0.03  um) 

y.mesh  loc=-0.84  spac=0.003 

#  Emitter  (0.05  um) 
y.mesh  loc=-0.79  spac=0.003 

#  Base  (0.55  um) 
y.mesh  loc=-0.5  spac=0.1 
y.mesh  loc=-0.24  spac=0.003 

#  BSE  (0.03  um) 

y.mesh  loc=-0.21  spac=0.003 

#  Buffer  (0.03  um) 

y.mesh  loc=-0.18  spac=0.002 

#  Tunnel  emitter  (0.015  um) 
y.mesh  loc=-0.165  spac=0.002 

#  Tunnel  base  (0.015  um) 
y.mesh  loc=-0.15  spac=0.001 

#  Window  (0.05  um) 
y.mesh  loc=-0.1  spac=0.01 

#  Emitter  (0.1  um) 
y.mesh  loc=0  spac=0.01 

#  Base  (3  um) 

y.mesh  loc=1.5  spac=0 . 3 
y.mesh  loc=3  spac=0.01 

#  BSE  (0.1  um) 

y.mesh  loc=3.1  spac=0.01 

#  Buffer  (0.3  um) 
y.mesh  loc=3.4  spac=0.05 

#  Substrate  (300  um) 
y.mesh  loc=303.4  spac=50 


b.  Regions 


#  Window  AllnP  (=InAsP) 

region  num=l  material=InAsP  x.min=-250 
0.84 

#  Emitter 

region  num=2  material=InGaP  x.min=-250 
0.79 

#  Base 

region  num=3  material=InGaP  x.min=-250 
0.24 

#  BSE 

region  num=4  material=InGaP  x.min=-250 

0.21 

#  Buffer  AllnP  (=InAsP) 

region  num=5  material=InAsP  x.min=-250 
0.18 


X .max=250 


X .max=250 


X .max=250 


X .max=250 


X .max=250 


y.min=-0.87  y.max= 


y.min=-0.84  y.max= 


y.min=-0.79  y.max= 


y.min=-0.24  y.max= 


y.min=-0.21  y.max= 


#  Tunnel  emitter 
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region  num=6  material=InGaP  x.min=-250  x.max=250  y.min=-0.18 
y .max=-0 . 165 

#  Tunnel  base 

region  num=7  material=InGaP  x.min=-250  x.max=250  y .min=-0 . 165 
y .max=-0 . 15 

#  Window  AllnP  (=InAsP) 

region  num=8  material=InAsP  x.min=-250  x.max=250  y.min=-0.15  y.max=-0.1 

#  Emitter 

region  num=9  material=GaAs  x.min=-250  x.max=250  y.min=-0.1  y.max=0 

#  Base 

region  num=10  material=GaAs  x.min=-250  x.max=250  y.min=0  y.max=3 

#  BSE 

region  num=ll  material=InGaP  x.min=-250  x.max=250  y.min=3  y.max=3.1 

#  Buffer 

region  num=12  material=GaAs  x.min=-250  x.max=250  y.min=3.1  y.max=3.4 

#  Substrate 

region  num=13  material=GaAs  x.min=-250  x.max=250  y.min=3.4  y.max=303.4 

c.  Electrodes 

electrode  name=cathode  x.min=-250  x.max=250  y.min=-0.87  y.max=-0.87 
electrode  name=anode  x.min=-250  x.max=250  y.min=303.4  y.max=303.4 

d.  Doping 

#  Window 

doping  uniform  region=l  n.type  conc=1.95el8 

#  Emitter 

doping  uniform  region=2  n.type  conc=2el8 

#  Base 

doping  uniform  region=3  p.type  conc=1.5el7 

#  BSE 

doping  uniform  region=4  p.type  conc=2el8 

#  Buffer 

doping  uniform  region=5  p.type  conc=0.95el8 

#  Tunnel 

doping  uniform  region=6  p.type  conc=8el8 

#  Tunnel 

doping  uniform  region=7  n.type  conc=lel9 

#  Window 

doping  uniform  region=8  n.type  conc=lel9 

#  Emitter 

doping  uniform  region=9  n.type  conc=2el8 

#  Base 

doping  uniform  region=10  p.type  conc=lel7 

#  BSE 

doping  uniform  region=ll  p.type  conc=2el8 

#  Buffer 

doping  uniform  region=12  p.type  conc=7el8 

#  Substrate 

doping  uniform  region=13  p.type  conc=lel9 

e.  Solving 
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method  newton  maxtraps=10  itlimit=100  climit=l  dvmax=0.05 
solve  initial 
solve  bl=0 . 9 

log  outf=multiIV. log 

solve  vanode=0.0  name=anode  vstep=0.025  vfinal=2.7 
extract  init  inf ile="multiIV. log" 

extract  name="Jsc"  y.val  from  curve (v . "anode" ,  i. "anode") 
X . val=0 . 0 

extract  name="Voc"  x.val  from  curve (v . "anode" ,  i. "anode") 
y . val=0 . 0 


where 

where 
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